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A METHOD TO ESTIMATE WEIGHT AND DIMENSIONS 
OF LARGE AND SMALL GAS TURBINE ENGINES 

E. Onat and G. W. Klees 



SUMMARY 


A method has been developed to estimate engine weight and major envelope 
dimensions of large axial flow aircraft jet engines and small gas turbine 
engines. The computerized method, called WATE-2 (Weight Analysis of 
Turbine Engines), determines the weight of each major component in the 
engine, such as compressors, burners, turbines and frames. A preliminary 
desian approach is used where the stress level, maximum temperature, 
material, geometry, stage loading, hub-tip ratio, and shaft mechanical 
overspeed are used to determine the component weight. 

A relatively high level of detail was found to be necessary in order to 
obtain a total engine weight within the required +10* accuracy. 
Component weiqht data for 29 different engines were used as a data base 
to develop the method for axial flow aircraft engines. The list of 
engines Includes military and commercial, turbofans and turbojets, 
augmented and dry, hardware engines and proposed engines, and supersonic 
and subsonic engines. For small gas turbines, weight correlations were 
provided through a subcontractor, the Garrett Corporation. Data used in 
those correlations were proprietary to Garrett and only the selected 
weight method is shown herein. The smal l-er,g1ne weight method is 
probably valid for other manufacturers, however this has not been 
verified. 

The accuracy of the method is generally better than +10X , on the order of 
+^*. The accuracy was verified by applying the method to 9 different 
engines, some of which were In the original data base and one small gas 
turbine engine. Engines used in the -validation were selected by NASA 
after completion of the program. 
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INTRODUCTION 


Aircraft anil propulsion system studies are frequently conducted by 
industry and Government. These studies may encompass a wide variety of 
enqine concepts ranqinq from relatively simple turhofans and turbojets to 
complicated variable-cycle enqines. The industry in qeneral has acquired 
an adequate computer capability to evaluate the thermodynamic performance 
of these diverse enqine concepts, however, an accurate method of 
estimation enqine weiqht and dimensions has not previously been 
available. The enqine manufacturers have developed suitable methods, 
however they are not available to the public. 

One method, 0^ that has been available to the qeneral industry 

predicted enqine weioht bv statistical correlations of major cycle 
character i st ics such as airflow, bypass ratio, overall pressure ratio, 
etc. This method is probably capable of rouqh estimates for conventional 
enqines; however, it is not applicable to nonermvent ional enqines and 
could not predict weiqht within + r > to *10t as would be required in 
typical preliminary desiqn studies. 

This proqram development was initiated to provide a more flexible and 
more accurate method based on correlations of component weiqht and 

physical characterist ics, such as compressor airflow size, pressure 
ratio, hub-tip ratio, etc. This type of approach then would be more 

capable of estimatinq nonconvent ional enqines, since the weiqht of each 

individual component would he accounted for. As shown in the Methods of 
Analysis, no adequate correl at ions could be found and a final nv.hod was 
chosen that is based on a mechanical preliminary desiqn which is 
responsive to the ma.ier enqine desiqn variables. 
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LIST OF SYMBOLS 



a 

A 

AR 

r 

Vyr 

C-0 

c/s 

0 

F 

g 

h 

h/t 

H 

HP 

J 

K 

L 

LP 

M 

N 

P 

PTO 

R 

RPM 

S 

T 

t. 

TLP 

TR 

U 

V 

U 


acceleration 

area 

blade aspect ratio 
Made chord, or convergent nozzle 
convergent-divergent nozzle 
solidity, blade chord to spacing 
ratio 
diameter 
force 

gravitational constant 
height, or specific enthalpy 
hub-tip radius ratio 
Total enthalpy 
high-pressure spool 
778 ft-lb/BTU 
factor for blade volume 
length 

low pressure spool 

mass flow, or Mach number 

number of elements 

pressure 

power takeoff 

radius, or gas constant 

revolutions per minute 

surface area, or blade spacing 

temperature 

thickness 

turbine loading parameter 
bladp taper ratio 
tip speed 
volume 

Weight, or weight flow rate 



t shear stress 

A turbine loading parameter 

0 Density 

u> Rotational Velocity 

0 Stress 

K ratio of specific heats 

e heat exchanger effectiveness 

0 ratio of local total temp, 

to std. temp. 

6 ratio of local total pressure 

to std. pressure 


Subscripts: 
h hub 

t tip 

b blade 

c case, or corrected conditions 

S stator 

hw hardware 

Stg Stage 

0 stagnation conditions 

bp blade pull 

? Engine inlet station 
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Attempts were made to improve this correlation. The data base 
information was examined to find if there were any mechanical design 
features (such as h/t, stage loading, etc.) that might account for the 
large weight difference. It was concluded that there were several 
different design characteristics that caused the lack of correlation, 
such as material properties, blade geometry, stage loading, shaft speed, 
and design 1 ife. 

A similar problem occurred for turbines (Figure ?), burners and diffusers 
(Figure T), augmentors (Figure 41, and duct heaters (Figure 5). Some 
components, however, did correlate well bv this method and these are 
further discussed in the Methods of Analysis. At this point it was found 
that the correlation method failed to predict many of the major component 
weights within the desired accuracy, and another method was adopted. 

The adopted method was based on a preliminary mechanical design approach 
where the design variables are taken into account. In the compressor, 
for example, rotor blade weight is calculated as a function of specified 
geometric parameters. Blade pull stress is then found, and the disc 
weight that will support the blade rotational force is determined. This 
type of preliminary design approach was also used for the other 
components (Fiaures 2 through 5). These methods are also discussed in 
greater detail in the following secton. 

The WATE-? method is intended to estimate weight for a given engine 
desiqn. It will not design an engine. This function must be performed 
external to the program. VIATE-? utilizes component state conditions, 
flow, work, etc., which are generated in the engine cycle portion of thp 
program (NNEPl. NNEP operation is described in reference 2 . 

In normal use of WATE-?, the desired engine cycle is simulated in NNEP at 
sea level static conditions for the engine design Doint. The user of 
WATE-? must be coqnizant of other conditions in the fliqht envelope of 
operation where maximum component temperature, work, speed, or flow 
occur. If these conditions are greater than the design values, they will 
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Figure 2. Data Bate Engines: Turbine Weight 
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sizo the component and will have a significant impact on weight. WATE-2 
allows input of scalers to account for these off-design conditions. 

A more accurate weight estimate can he achieved by developing an array of 
engine cycle data over the intended flight envelope. The WATE-2 program 
will scan the output data of the engine cycle and select the maximum 
conditions for each component. This automated procedure eliminates the 
need for scalars to determine maximum component conditions. 

The NNEP simulation of the engine may require the use of components 
are required mathematical ly, but are not desired in the engine weight or 
dimensions. These can be selectively eliminated. 

The methods of analysis described for each component in the following 
section have been developed to achieve an overall accuracy of ±10*. 
Since the rotating components comprise the major part of the total 
weight, considerable detail was necessary in order to achieve the 
accuracy goal. Normal program users may not have sufficient knowledge to 
adequately define all of the necessary inputs; however, typical values 
are given in the User's Guide, Section 2.0. 

1.1 Large Gas Turbine Engines 

The methods developed in this section were based on correlations of the 
data-base family of large gas turbine engines. Table 1. There are 
certain components that should be equally applicable to small engines 
such as: annular ducts, burners, rotating splitters, shafts, mixers, 

nozzles, certain frame types, transmissions or gear boxes, and heat 
exchangers. Most notable differences between large and small engines are 
the use of radial-flow components, and relatively thick discs for small 
axial-flow components. This section should otherwise apply equally well 
for both large and small engines. 
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1.1.1 Axial Flow Fans and Compressors 

The general procedure used for fan and compressor weight prediction is a 
stage-'oy-stage mechanical design as illustrated in Figure 6. Rotor blade 
volume and weight is determined; then, blade pull stress, disc stress and 
disc volume are calculated. Connecting hardware, stator, blades, and 
case are then estimated and summed to give the total component weight. 
The following input data are necessary: 

o An allowable pressure ratio for the first stage which reflects 
the design approach and technology level. Specific work for 
this stage will be held constant for additional stages. Number 
of stages can also be specified as an option. 

o The entrance and exit mach number of the component. 

o The hub-tip ratio of the first stage. 

o Compressor design mode: constant mean-1 ine, constant-hub, or 

constant-tip diameter. 

o Effective density of blade material: defined as total blade 

weight divided by total volume. 

o Maximum inlet and exit temperatures, if not at design. 

o Aspect ratios for the first and the last stage blades. 

0 N max/ N des overs P eed factor, 

o Blade solidity, 

o Density of disc material, 

o Blade taper ratio 

o Blade volume factor, ratio of total volume to blade volume 
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diameters of the flow path are then calculated from the specified radius 
ratio: 



4A 

ir(l-TT7 2) 


(1) 


0 


h 



( 2 ) 


Compressor RPM Is determined hy dividing the tip speed (found from 
Figure 7) by the stage radius (R t ). 

Dimensions of succeeding stages are based on the design mode selected 
(constant mean, tip, or hub design). In the constant mean-line method, 

the mean radius is based on the mean flow area of the first stage. 

Corrected airflow at the entry of each stage is determined from 

calculated state conditions. Stage entry Mach number is assumed to vary 
proportionally to the number of stages when a different entry and exit 
Mach number are specified. Tip speed for the downstream stages are then 
calculated from the stage dimensions and shaft speed. Radius of the 
mean-area location of other stages is set equal to the first stage 

mean-area radius. 


Blade aspect ratio Is determined by assuming a proportional change for 
each stage If entry and exit aspect ratio are different. Volume of metal 
in the blades is then calculated by: 





O' 



Where K Is a volume factor which accounts for fir-tree mount volume, 
taper ratio, and thlckness-to-chord variations In typical blades. For 
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the data base engine, K .as found to be 0.055 for fan blades, .1? for 
compressor blades and for blades with h/t greater than .8: 


K - .120 + .04 ( h/ t - .80) 


(4) 


The rotating blade might of each stage Is determined from the blade 

volume and materia, density specified. Haterla, density automatical,; 

changes from r, to steel .hen the stage temperature exceeds a specified 
maximum, normally 700°F . 

* . ’ " , 1371 CK sta « e temperature Is calculated 

from the NNEP output data and the specified over-temperature ratio. 

a or wight and dimensions are assumed to be equal to the rotor blades 
and include the inner diameter stator shroud. 

Stage length 1s found In the following manner. Blade chord (r ) as 
shown in Figure 8, Is the quotient of AR and the blade height: 8 


(5) 


r - n 
C B ' “A F“ 


I 


tor the data base engines, the stator length .as found to be egua, to the 
rotor length (or blade chord,. .1th ,7* of the rotor length required for 

s uLTnd T ^ SUt0r " d the between the 

7 ' "" S ” ° f *" “** s »«"• the ‘»t„ compressor 

c d ’I T !V 9 VW " ar<! " 0t ,nC,Uded ,n th ‘ compressor wight, but 

can be Included as a frame, see Section 1.1.7. 

Total number of blades Is found from the specified solidity (C/S) and the 
chord length determined above: 


N 0 « 2ir • R t • C/ ♦ AR 


(«) 
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This value is truncated to an Integer number of blades, and the same 
value is used for the stator. 


The maximum blade-pull stress occurs at the blade root and is a function 
of tip speed, blade height, and material density. Expressed In terms of 

the nondlmenslonal Input blade geometry, the equation for blade-pull 
stress Is: 


12PU, 


. h/ t 2 


* TT 0 - h/ t) (l + 3 h / t )] 


where U t is the product of tip radius, N and the overspeed factor 

Woes' 

The compressor discs are a larqe part of the total engine weight and It 
is therefore necessary to define this weight as accurately as possible. 
Oise weight is a function of: diameter, blade load carried on the rim of 
the disc, material type, speed, disc shape (or thickness distribution), 
and design stress level selected for life considerations. A theoretical 
functional correlation was developed which showed that disc thickness 
should vary inversely with the product of blade-pull stress and disc 
diameter, i.e., the blade load per unit thickness. For those engines in 
the data base where large-scale drawings were available, several typical 
discs were measured. Blade-pull stress and disc volume were calculated, 
and the results were plotted in the form of relative disc thickness 
versus load per unit thickness. Figure 9. 

There is an identifiable trend in these data that show a different 
characteristic for compressor discs and fan discs. Low hup-tip ratio of 
the fan probably accounts for the departure from the compressor trend. 
The allowable stress difference between steel and T1 causes the different 
trend for the two materials. Further Improvement* in the accuracy of the 
disc volume correlations could be made If the number of discs were 
increased In the data base and the rim loading term was modified hy the 
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design Stress ,eve, or maxim™ a„owab,e. An options, disc preHminary 

Z'* r '* aU0 avai,ab,e a " d recommended tor pote 
improved accuracy (see sec. 1.2.3). potential 

the^ata^ase UStrateS “ SU9e ”' < ' th0d * hat was used <" «<*t of 

; 5 r ; "sc d Cer “ tyPlCa "' V 3 •° 75 - ,n thiCk Cy "'" der ’taw at 
_ radius. The connecting hardware (u ) is estimated by the 

following equation: hw y 


W HW = .75 x 2. x irR n x .075 x L $TG x p 


(a) 


where R h is the disc radius (or blade hub radius) L is th* 
length, and p is the material density. ’ Stg 9 * 

i:;r^r I s the iast item ° f wei9ht inciuded in the compress ° r 

oht buildup. Average case thickness in the data base engines was 0 i 0 

calcXTd "t T""’ 1nClUdin9 ^ flangeS ‘ Ca$e 

assul Lr thTcaL: 9 " ^ ^ «* d ’“ * 


w c = * D t x l stg x - 10 x P 

where D t is the stage tip diameter. 


(9) 


ZLTZr i9ht 15 S “" ° f r0t0r b,ade - Statdrs - d ’«. connecting 

weight. ”* ^ * WG19htS *"* SUmmed t0 * ive the total component 

i ner tiiTof JacTst ^ 3 " d turb1nes is determined by finding 

each stage of each component, it is assumed that hlades have 

uniform weight/area and the disc is of uniform thickness. This method 
een 3PPlled t0 severa1 existing engines and the results showed good 
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agreement. The following equations are used for the blades and disc 
inertias: 




f 

! 

m 



a 


I 


WT BL h B 2 R 2 x N 
BLADE g x "1 + h B 


! DISC * ^ jp - SC x R h 2 


The total stage inertia is the sum of blade and disc inertias 


^TG " 1 BLADE + I DISC 


The component inertia is the sum of the Inertia of each stage. 


( 10 ) 

(ID 


(12) 


1.1.2 Axial Flow Turbines 

The method described for compressors is generally the sam for turbines. 
Input data required are: 

o Maximum tip diameter of the first stage, or number of stages 

o Inlet Mach number (axial) of the first stage, and exit Mach 

number (axial) of the last stage 

o Rotor blade aspect ratio of the first and last stages 

o Solidity 

o Reference disc stress, 0.2* yield point of the material selected 

o Cooling indicator, to modify the blade volume calculation for 
cooling holes 
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o Design mode, constant hub radius, constant mean area or constant 
tip radius 


o Shaft overspeed factor 


o 


o 


Turbine loading parameter. 
Blade material density 



o Blade taper ratio 


o Blade volume factor 


Two alternative procedures can be used to size the turbine: (1) specify 
maximum diameter of the first stage and find the number of stages from 
the work loading parameter, or (2) specify the number of stages and the 
diamter of the first stage is found from the work loading parameter. 
Shaft speed is transferred from the corresponding compressor; and in the 
case of (11, the number of stages is found by iteration until the 
resultant tip diameter is equal to or less than the specified diameter. 
Equal work per stage is assumed. Total component work and state 
conditions are taken from the NNEP stored cycle data. 

To determine blade height, the flow area necessary to pass the corrected 
airflow is calculated at the inlet of each stage. When the first stage 
inlet Mach number is different from the last stage exit Mach number, a 
proportionate change is assumed for inlet Mach number of the other 
stages. Hub radius of the first stage is found by subtracting from the 
stage-projected area the exit-flow area required to satisfy the specified 
exit Mach number. Dimensions of the remaining stages are then determined 
from the design mode specification and the calculated exit area. 



Aspect ratio and number of blades 
in determining the compressors. 
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same method. Equation (3), except that K * 0.195. When the blade is a 
cooled blade ( normally HP turbine blades with relatively low-aspect 
ratio), the calculated volume is reduced 20% to compensate for 
cooling-air passages. Blade weight is then found from the specified 

material density. Blade-pull stress is calculated by Equation (7). 

In a manner similar to the compressor discs, turbine discs were measured 
in the data base engines to produce the results shown in Figure 11. Th-' 
rim loading parameter (o Bp x R n ) was modified by dividing by the 0.2% 
yield strength (at the max. operating temperature) of the particular 
material used in each disc. There are many different steel alloys and 
superall oys that have different strength capabilities for which this 
correction compensates. Each disc has a different operating temperature, 
maximum allowable stress, and design life. These factors and other 
unknowns, such as the effect of hub tip ratio, are believed to be the 

primary cause of the data scatter. The best-fit curve shown in Figure 11 
is used in WATE-2, unless the optional procedure is used where the discs 
are calculated by a preliminary design procedure developed primarily for 
small engine discs. This method is recommended because it includes the 
effects of blade aspect ratio and disc stress level. 

The relative disc thickness is found from Figure 11 with the calculated 
independent variable (°gp x R n ) . Disc volume is found by multiplying 
the relative thickness parameter (V/D h 2 ) by D h 2 . Blade-material 
density is an input; however, disc material is assui id to be steel or 
superalloy with 0.286 density. Since all of the data base engines used 
steel or superalloy discs, correlations of other materials could not be 
j made. 

Each stage of the turbine is treated as a stator-rotor pair (as opposed 
to rotor-stator oair in the compressor. Stator blades are also assumed 
to have the same rurnher, volume of material, and dimensions as the rotor 
blades. The stator weight is calculated ^y Eauation (3), with K=0.144. 

Stator-rotor spacing is the same as compressors, 17% of the rotor 

length. Stator AR is taken to be 83% of the rotor blade AR. 
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Connecting hardware and case weight are also determined by the same 
method as used In the compressors. The total weight and length of the 
turbine component Is the sum of disc, blade, stator, connecting hardware, 
and case. No exit guide vanes (E6V) are assumed In the turbine 
component. EGV's, If required, can be considered a part of the exit 

frame weight. Rotational Inertia Is determined In the same manner as the 
compressor. 

1.1.3 Ducts 

It is assumed that the major structural load In a duct Is a result of the 
Internal pressure. Also, It Is assumed that the Inner wall of the duct 
is the same gage as the outer wall. In the data base engines the outer 
surface of the OD walls were typically exposed to ambient pressure. The 
ID wall was subjected to fan pressure, hp compressor exit pressure, etc., 

and aP for the ID wall could not be generalized. Figure 1? Illustrates 
the duct and nomenclature. 

The equation for stress ui a longitudinal section of a thin-walled 
cylinder subjected to an Irternal pressure is : f 3 ' 

p °n 

0 * (13) 

and for solving for minimum thickness: 

P D q 

*ta1n = To (14) 

where o Is the allowable for the material, p Is the Internal total 

pressure, and Dq is the duct outside diameter. 

T1 ^is assumed with SO, 000 lb/1n P allowable at temperatures below 
700 F, and steel Is assumed at 70,000 lbs/1n ? above this 

temperature. The appropriate material Is selected based on the total 
temperature of the duct airflow. The weight Is calculated as a function 
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duct length (L) f the Inner diameter (0.), and the outer diameter 

<V : 

“tlUCT -’'< 0 o + U 1> LP ‘»1l» (15) 

Duct Mach number Is specified as an Input, and corrected airflow Is 
determined from the NNEP cycle data. ID and 00 are determined as a 

function of required flow area and the dimensions of the connecting 
upstream component. Care should be taken to ascertain whether these 
assumptions apply for specific engine configurations. For example, a 
thin-walled cylinder subjected to an external collapsing pressure will 
fall at a much lower pressure than It would If It were subjected to an 
Internal bursting pressure, as assumed In the duct weight calculation. 
If both ID and 0D walls of the duct are exposed to ambient pressure, the 
10 wall should be sized to avoid collapse, such as determined 
experimentally by Stewart 3 for lap-welded steel tubes: 


P MAX = 1000 [ 1 * V^-ISOO (§) 2 ] (16) 


or expressed in terms of minlmtan wall thickness: 



(17) 


MATE-? does not perform the above calculation to determine whether 
collapsing pressure sizes the 10 wall. It also does not check to 
determine whether the gage selected Is less than minimum gage. 

1.1.4 Rotating Splitter 

A rotating splitter, Figure 13, Is a circumferential separator of two 
flows within the same compressor. These flows normally have different 
pressures and temperatures, and the splitter must perform a sealing 
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funaion. Stages that Incorporate rotating splitters are treated the 

^nnrr : ,'rrrt: , i — - «-• — 

rotor-stator spacing Is 17* of rotor length." 16 ^ b ' ade ' 

centrifugal "blade' f‘ ter T *° ** b,a<,e and the 

the added load. Sp.Tuer ,£*, T ^ 

ofiuiier weight, per blade Is estimated by: 


C ’ x .10 


— » »> - >> - ~ 

rlr " ” ™ ~ *>•"'« « 

r:— r;x: - 

The centrifugal force contribution of the rotating splitter is: 


F - MrW 2 = x R x ( & x RPM \ 2 

9 s V 6 


:• * - — 

p stress (Equation 7) is Increased by the amount 


AO . ^PL. 
BP tf- 


x t/c 


where C fi Is determined by Equation («5) and t/c is th* thin. 
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using the disc 


volume 


correlation. 
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urner flow area is determined from the input velocity th* m - 4 - 

and the entry-corrected airflow from the NN£P cycle data ^ ’ US ’ 
Is found to give the sn»rifi*n cyc,e data - Bur, » r length 
and entry renditions. time based on the Input velocity 


1 ‘ ,,n im6r 0^ the bd™r 

„t a " d V “Ith the specified mean radius located at 1 
Outer-case thickness Is determined by Equation iU) Th „ " 

Is used for the inno. equation i!4). The same thickness 

io.ooo r ™ 6 ,s •*-> -«■ 

inner and outer cases are fn, A c ' VOlunie ° f mater1al f or the 

n V Equat1on ( 15 ) • using burner length (L). 

st"o e i wTnt: LtTZ'Z a me manner - 0 - 055 - ( " *"ck 

case. The burner dome, fuel ,7 othT a " d ° Uter 

are estimated hv tvolcal oonma f * , * d th com Ponents 

determined by the following equation^ *" ^ ^ ^ ^ * S 


DOME 


■ .0106 (R t 2 - R } 2 \ 


( 20 ) 
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Total burner weight is the sum of inner and outer cases and liners, 
burner dome and fuel nozzle system, and frame where applicable. 


1.1.6 Shafts 

A shaft is assumed to be the power connection between components, see 
Figyre IF. Multiple stages within a compressor or turbine are also 
connected by a shaft; however, this weight is included in that 
component's weight buildup. 


The required inputs are: 

o The component numbers connected (to determine length, power 
transmitted and shaft speed). 

o The shaft material density and allowable stress. 

o Radius ratio (of the inner shaft only). 


Multiple concentric shafts can also be specified, and will be sized 
around the inner shaft with 0.20 in radial clearance assumed. 


Dimensions of the inner shaft are determined to provided the necessary 
torque capability at the specified allowable stress. Total shaft power 
is the summation of work (A Hi for all turbines on the shaft. Torque is 
calculated by: 


T - — x J 

fa) 


( 21 ) 


where u is the shaft ;otative speed, rad/sec. 


Shear stress due to the torque load is defined by ( 3 ) t 


i 


16 T D q 

’ » <°„ 4 - » 1 4 > 


( 22 ) 
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or In terms of the Input diameter or radius ratio (°^/D 0 ) : 


16T 


Solving for D Q in terms of allowable stress (o ) : 

a I ! 


D_ = 


Weight is then found by 


16T 




1/3 


(23) 


(24) 


W * L p ir D o 2 [ 1 - 1 

T~ L D. 2 J 


(25) 


A similar procedure is used for concentric shafts. The second shaft's 

inside diameter is found by adding 0.40 in to 0 Q , and Equation (24) is 

solved by iteration. 

While it is assumed in the shaft-weight estimate that torque determines 
the shaft dimensions, it should be recognized that other design 
considerations may dictate shaft dimensions. Shaft critical speeds or 
longitudinal stiffness may actually design the shaft, but this is a 

function of bearing arrangement, mount stiffness, location of and 

stiffness of rotating masses. The calculated shaft weight should be 

considered to be an absolute minimum, and can possibly be much larger 
when these other criteria are considered. The rotational inertia of the 

shaft « not calculated, since it is a negligible quantity compared to 
the compressors and turbines. 
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1.1.7 F rames 

A structural frame is formally required to span the enqine flow path from 
the outer enqine case to the shaft, usually to support a bearing (as 
shown in Figure 16) for several typical engines. Mechanicel design of 
the frame would require a definition of all loads imposed on the frame 
under normal operating conditions, transients, and adverse operating 
conditions, such as a hard landing. This level of detail is normally not 
available at the preliminary design stage for which WATE-2 has been 
developed. 

It has been found, however, that the frame weight of the data base 
engines correlates well with the total frame-projected area. These data 
are shown in Figure 17 for four types of frames commonly used: 
single-bearing frames with and without PTO, turbine exit, and 
intermediate. Frame weight is determined from these data, based on the 
local diameter and the type of frame specified. 


1.1.8 Nozzles 

Unlike the rotating components, the loads and load paths of nozzles 
(particularly variable area C-D nozzles) are not readily defined on a 
general basis. A selected type of nozzle could be subjected to a 
detailed weight-estimating procedure, however, the trade-offs of internal 
and external performance with nozzle length and diameter would also be 
necessary to optimize the design. This type of data is not likely to be 
available at the level of development for which WATE-? is intended. 

A procedure has been developed that shows proper trends for 

multiple-stream nozzles and for variable geometry and fixed-geometry 
nozzles. Nozzle length is specified and should be selected to be 
representative for the type of nozzle, i.e., C or C-D. An effective 
surface area is calculated based on the diameter of the connecting 
component and the specified le -th. 
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Figure 17. Frame Weight 








Only circular, conical nozzles are assumed, however, coannular nozzles 
could be represented by specifying a circular nozzle for each flow path. 
Plug nozzles can be represented by specifying a larger effective length; 
e.g., from nozzle entry to end of plug. Wall thickness is assumed to be 
0.10-in steel above 700°F and Ti below 700°F. Variable nozzles are 
calculated in the same manner except that the effective wall thickness is 
2.75 times that of the fixed nozzle. 

1.1.9 Mixers 

A mixer is a device placed at the point of confluency of two coannular 
streams to increase the mixing boundary so that thermal mixing takes 
place in a minimum length. This type of mixer is sometimes called a 
daisy-mixer, chute-mixer, or forced-mixer. 

Flow area of each annular path is taken from the NNEP cycle data, and the 
inlet radius (Ri) of the upstream component is used as a starting point 
for locating and Rq, as shown in Figure 18. Normally Ri will be 

the hub radius of the last turbine stage. 

Mixer length (L) and number of passages (N) are required inputs. The 
following relationship has been developed that is representative of the 
surface area of typical mixers: 

S = [3.93 R m + 1.25 (R Q - R^] L (26) 

material assumed is 0.10-in thick steel. 


1.1.10 Annulus Inverting Valve (AIVl 

This device has been used in some variable-cycle engines to invert the 
annular position of two concentric flow paths. It accomplishes the flow 
inversion within a constant diameter envelope, and with constant-area 
duct passages. Figure 19 shows a typical example of an AIV. This AIV 
was designed to vary the byoass ratio in a JT8D engine. 
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RM IS BASED ON INNER AREA AND R, 

Lgp = L 2A tw ~ INPUT 
A • AquTER * A INNER 
L IS DETERMINED 

WTM - (K, * R M ♦ 126 x N (RO - Rl»x LxKj 
K,- 3.927 Kj * 028 

Figura 18. Mixar Schematic 


43 


* 



The AIV weight method assumes e coostruction similar to that shown In 
-gore ... except that instead or sheet-metal, the materia, is assumed to 

87 Jn'T “o' ,h/,t be '“' 700 ° F »"« honeycomb at 

1.87 lb/ft above 700°F. 

specified. 


' "VMVJV.U1IIU 

If desired, different materials can 


AO empirical relationship, similar to the mixer method, has been 
developed for estimation of the AIV weight: 


U 


AIV 


= L x 


}f [* » R, + 2w R 0 + 3.93R^ t 1.25 (R^)] 


where R, ,s the hob radios of the opstream connecting component, and 

* a " d R 0 are found t0 satisfy the inpot Mach nomber with the entry 
corrected airflow (see Figure 20). The number of passages (N) is an 
input, and material weight per square foot (N/A) is selected depending on 

the stream temperature. Length (L) of the AIV is calculated from the 
input specific length, L sp : 

L 


-SP 


\/4A/n 


(28) 


Specific length is preferred as an input because it is nondtmensional. 
and it is a major variable that determines AIV pressure loss A 
relatively good compromise between sire and performance is achieved when 


N>-1 and L SP =0.R to 1.0, which results in a pressure loss between 7 5* 
ind 1.M5. 


If the AIV is of the switching type, where one half indexes in a 
rotational direction relative to the other half to change flow-path 
orientation, an actuator weight is estimated at 10X of total AIV weight. 
Additional structure to support the rotating half is not incuded an d 
should he represented as an additional frame. 
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Lgp - L/ViA/ir - INPUT A - A INNEr| 

FROM MACH INPUT 

a INNER AND ^UTER ARE found 
L IS DETERMINED 

WTIC-2 R, L WTSI 
WTOC * 2 • Rq L WTSO 

WTWALL - (K, - R M + K 2 N (R 0 - R, )) L WTSW 
WTSI. WTSO AND WTSW ARE MATERIAL WEIGHT IN lb/ft 2 
R M - ^A/IT+R, 2 

K, -3.927 K 2 -1.25 


Figure 20. Annulus- Inverting Valve Schematic 
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*11 Transmission 


A method of estimating the weight of various types of gear systems has 
been previously developed by Schmidt. 4 This method is used in WATE-2 
for gear sets that are typically used in aircraft. 


The equation for the total gear box weight is: 

“ ■ * 4 ( 3 *) - 80 


(29) 


1.1.1? Thrust Reversers 

A weight estimating method previously developed for aircraft preliminary 
desiqn studies is based on the weight of 18 different reversers that are 
in current use. It has been found in correlations of these data hat 
reverser weiqht (W) is a function of corrected mass flow anc j 

nozzle oressure ratio fP R l, and is dependent on whether the iream is 

hot (primary) or cold (fan). The following relationship has been 
developed: 

[K,XW^ + K 2 ][K 3 P R fK 4 ] f30) 


where hot streams K, = .52. K ? = 423. K . 1.004 and K. . -.5054 
For cold streams K, = 2.22, K ? = 11.0, K, . .23, and . .56. 

The MATE -2 method will apply the cold stream equation to a fan stream 
whether or not it is heated by a duct burner. The hot stream equation is 
used for turbine outlet streams or mixed-flow exhaust streams. 

1 . 1 . Heat Exchang ers 

Both rotary and fixed heat exchangers can be estimated. Methods 

previously developed produce adequate results for preliminary design 
purposes, see Figure ?1 . 
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For rotary heat exchangers, a ceramic core is assumed. Weights of this 
type ofheat exchanger have been determined by the Corning Glass 
Company* and are represented in Table 2 for various levels of 
effectiveness and pressure loss. These data are developed for a total 
corrected airflow of 200 lbs/sec. For other sizes, * weights are 

scaled directly with corrected flow. 

Fixed-tube heat exchangers are estimated by ? eat transfer analysis* 
whereby the required tube surface area is f . d to give the specified 
effectiveness. Flow area of the tubes i* founv' from an input Mach 
number, number of tubes, and corrected flow. 

Wall thickness of the tubes is determined by Eouation (141 to satisfy an 
assumed allowable stress of 50,000 lh/in? and a density of 0.168 below 
700°F. A stress of 70,000 lb/in 2 and a density of 0.286 is assumed 

above 700°F. The length of tubes is determined to satisfy the surface 
area requirement. 


7a/)/g 2. Ceramic Rotary Regenerator^ Weight 



i 

200 lb/$«c corrected wtiyht flow 

BPR 

3 

7 

10 

Ap/p (%) 

5 

10 

5 

10 

5 

10 

Weight (lb) 







e • 80% 

674 

642 

628 

446 

600 


86% 

971 

754 

914 

640 

834 


90% 

1,622 

1,188 

1.428 

1.085 

1,313 

994 
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Fixed- tube heat exchanger tube weight (W t ) is then found b.v 

W t* PL t w (R 0 2 " R 1 2) (31) 

where R and R- are the tube radii and L is the total length 
oi . 

required. Casings, mounting hardware, manifolds, and other equipment 

that may be necessary are assumed to be equal to W^. Total heat 

exchanger weight is 2W^.. 

\.2 Small Gas Turbine Engines 

Discussed in this section are all of the components that are unique to 
small gas turbine engines. These include centrifugal compressors, radial 
turbines, axial compressors, axial turbines, burners, cross-over ducts, 
diffusers, compressor shroud, aear box, and rear frame. Some components 
utilize the same methods as used for the large engines, and these are 
noted in the discussion. A typical small-enqine is shown in Figure 
where the major components are identified. 

i . ? . 1 C entrifugal Compressors 

The following input data are required for centrifugal (radial) 

compressors: 

o Entry axial Mach number 

o Max first st.aqe pressure ratio 

o Radius ratio of the inlet hub to the inlet tip 
o Shaft overspeed factor 

If the specified overall pressure ratio is greater than the max allowable 
first stage pressure ratio, another stage is added and the work is 
divided equally. Stages are added until the first stage pressure ratio 
is loss than the max allowed. Tip speed of the second stage of multiple 
stages is taken to he QpM of the first stage tip speed. 
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gearbox 


Figure 22. Small Engine Components 
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Hub and tip radii of the entry are found from the input radius ratio, the 
entry corrected airflow, and the entry Mach no. The tip radius of the 
outlet of centrifuqal compressors is assumed to be 3 times the hub radius 
of the entry, which is typical for wheels from 10 to 20 in. diameter. 

Corrected tip speed is found from the calculated pressure ratio of the 
staqe, Fiqure 23, and shaft speed is then found from tip radius and tip 
speed. 


The shaft speed will normally be set by turbine stress. Since turbine 
stress is a fall-out calculation in WATE-2 (for axial turbines only), it 
should be checked to see if it is reasonable, and the input radius ratio 
should be adjusted if necessary to chanqe the calculated shaft speed. 

Centrifugal compressors are assumed to be titanium and- the weight is 
found as a function of tip radius, Figure 24. Length and rotational 
inertia are found similarly, Figures 25 and 28. Figure 28 also compares 
the inertia of radial turbines. 

1.?.? Radial In-Flow Turbines 

The tio radius of radial turbine wheels is found by the inlet corrected 
airflow, Fiqure 27. The technoloqy level, current or advanced, must be 
specified. Length of the impeller is found from Figure 28, as a function 
of the tip radius. 

Steel turbine wheels are assumed, and the weight and rotating inertia are 
determined from t i d radius, Fiaures 29 and ? 6, respectively. 

1.2.3 Axial Compressors and Turbines 

This method is identical to the larqe gas turbines except for the disc 
weiqht c •' dilation procedure. No additional inputs are required, except 
for the input which activates the optional disc procedure. It was found 
that the previous disc weiqht correlations were not applicable to small 
diameter discs. Typically, the small discs were 3 to 4 times heavier (or 
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Figure 24. Weight of Titanium Centrifugal 
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Figun 26. Rotating Inertia of Centrifugal Compreuon and Radial Turbinm 








thicker) than would be predicted by the original method. A new method 
was developed specifically for small -engines. It is also applicable to 
large engines, and may provide improved accuracy. 

An inserted blade is assumed, and the disc material is assumed to be Ti 
or steel forgings with ultimate strengths of 1?0,000 and 160.000 psi 
respectively. The rim of the disc is assumed to be 10* of the blade huh 
radius or .75 in, whichever is greater. The remaining disc volume is 
assumed to be a trapezoidal section rotated about the axis. T e 
tbi«-kness of the trapezoid at its outer radius is sized for 75* of t e 
ultimate strength. Thickness of the trapizod at the inner radius is 
bas*d on 50* of ultimate strength for the tangential stress. This va ue 
„ selected based on an experimentally determined hurst speed marptn. 

,,,0 average tangential stress can be satisfied with a constant 
thickness disc, the here radius is increased until the lowest weioht disc 
is achieved that .iust meets the design criteria. 

stresses are calculated based on the total blade force act in, on the 
outer circumference nf the disc, and it is therefore sensitive to chances 
in made aspect ratio and solidity as well as blade pull stress and huh 
radius. The original method is not sensitive to aspect ratio or solidity. 

1 - .1 Crossover nuc_t_ jnt eqr a 1 with _P iffuser 

IMS component TVTnlv used to connect two centrifuoal compressors and 
„ has an inteural diffuser. The diffuser outside diameter (0 is 

h( , 1.S tunes the first impeller diameter and the length (1 is 
he of the diff»«~ di, meter. *>«** of the crossover duct 

is osMmit.’d Kv tho following relationship: 


VJT - 11.5 L x 0 


(T?) 


o 


whi’fo l and n am in ft 
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1.2.5 Diffuser 

A diffuser is required when only one centrifugal stage is used. It is 
normally feeding into a burner directly from the compressor. The outside 
diameter of the diffuser is also 1.8 times the tip diameter of the 
impeller, and the length is 7X of the diffuser diameter. The weight is 

expressed by: 


WT = 32.8 L x D. 


(33) 


1.?.6 Rear Frame 

For small gas turbines, a turbine exit frame is normally required to 
supoort the turbine bearing, radial loads only. The front frame and 
front bearing take the thrust loads, and this is sometimes an integral 
part of the gearbox. The following equation is used to estimate the 
weight of the rear frame, which includes the tailpipe and nose cone: 


WT » 55.5 R t 2 + 6.53 


(34) 


where R t is the turbine tip radius (ft). 

1.7.7 Centrifugal Compressor Housing 

This item is similar to the outer case in an axial compressor. Weight is 
based on the length and the tip radius (R t ) of the impeller: 


WT = 163 x L x R t 


(35) 


1 o a Gearbox and Accessories 

The existing transmission weight method will give satisfactory results 
for small enaines. Accessory weight however is a significantly higher 
percentage of total weight. Accessories are typically about 17* of the 
total engine weight, excluding accessories. 
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Other Program Function s and Capabilitie s 
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maximized values are s scanned over all 

The ™;r nigbrXl and" the^maximum conditions are 

.. *«• «• - — - - the 
maximum condition occurs for each component. 

.rP^fth. user to achieve a reasonable engine design, the 
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r:.ir. .... . 

shows a list of warnings and corrective actions. 
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The WATE-? proqram will automa ic elected scale 

. c • si in thp NNEP simulation. Up to . 
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table 3 design limits 


WARNING MESSAGE 

ACTION RECOMMENDED 

r 

OEFAULT VALUE 
TESTED AGAINST 

BLAOE Pull STRESS EXCEEDED 

REDUCE SHAFT SPEED 
(BY RPM SCALAR) OR 
INCREASE EXIT MACH 
NUMBER 

O « 50000 HPT 

0 = 60000 LPT 

0 * 80000 FAN & HPC 

h/t TOO LARGE 

REDUCE HUB-TIP RATIO 
INPUT 

h/t = .93 HPC EXIT 

h/t TOO SMALL 

INCREASE HUB-TIP 
RATIO INPUT 

h/t = .32 FAN, HPC ENTRY 
8 .50 HPT, LPT EXIT 

STAGE LOADING TOO HIGH 

INCREASE A INPUT 
(WHICH IS 1/LOADING) 

A = .28 HPT, LPT 

STAGE PR TOO HIGH 

REDUCE STAGE PRES- 
SURE RATIO INPUT 

PR = 1.8 FAN 
8 1.4 HPC 

FLOW VELOCITY TO HIGH 

DECREASE STAGE INLET 
MACH NUMBER INPUT 

M ex = - 60 ALL COMPONENTS 

BLAOE SIZE TOO SMALL 

CHANGE DES OVERALL 
PRESSURE RATIO OR 
REDUCE H/T INPUT 

hg .40 ALL COMPONENTS 




relative to the front-flange of the engine is determined as a function of 
the position of each component. 




For rotating components, the weight of each stage is assumed to act at 
mid-length of the stage. Moments are summed about the front flange for 
each stage. The center of gravity of each rotating component and the 
total engine is calculated. CG locations are shown in the output data 
for the total engine and each component. 

1 .4 Program Val i dati on 

A verification of the accuracy of the method can only be done by applying 
it to various tvpes of engines and comparing the results with the actual 
measured engine weight and dimensions or with those estimated by the 
manufacturer for Droposed engines. Since the manufacturer's estimate of 
proposed engines also includes some error, the real deviation or error of 
the WATE-2 method can only be found by comparing engines that have been 
built in production quantities. 

In order to judqe the accuracy of the method, the NASA program manager 
se -| e cted 3 large engines and one small engine for comparison. These 
included both production and proposed engines. The selection was made 
after the method was completed and submitted for approval. Results of 
the WATE-? estimates for these engines are shown in Figure 30. As can be 
seen, dimensions and weiqht of the 0 selected engines are within the +10* 

accuracv goal . 
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2.0 USER'S MANUAL 

Ihis section contains a description of input-output data, values of 

VP1 1nDUt$ 3nd $ample cases * WAT E-2 is designed to function around a 
component type engine cycle analysis program, the "NASA-NAVY Engine 

ogram. Reference 7. The calculated thermodynamic output data of NNEP 

not described here since it is unchanged from its normal mode of 
operation. 

The thermo c .sign point case of NNEP can be used to generate the WATE-? 
.nouts or additional NNEP off-design points can be run and the output 
ata -HI be scanned for maximum conditions of shaft overspeed, work 
ow or temperature for each component. In order to produce the most 
accurate weight estimate, the pff-design cases should encompass the 
maximum performance level required of each component. All components 

! : 0nt ;: bUte we ' 9ht ™ St a ' S0 be inc ' ude<l the NNEP enqine model. 
WATE-P will not calculate weight for components which are not included in 
the NNEP engine simulation. 

WATE-P also Will accept an input weight scaler for each component so that 
selected components can be increased or decreased for eliminated 

altogether! to determine sensitivities, etc. Both SI and English units 
of measure can be used. 


2.1 Program Structure and Deck Stacking 

The overall prooram structure is shown in the flow chart. Figure PI . 
NNEP design point data is stored In the "thermodjnamic data" for use in 
calculation weights and dimensions (W/Dl in the WATE-? part of the 
program. Computer execution time, core storage, and cutout print 
requirements have been increased s'iqhtlv over the NNEP program. 

The order nf deck stackinq of the ,lob Control cards and NNEP input data 
are unchanged from the normal operation „f NNEP. Two new inputs have 
been added to the NNEP indicator set. The indicators, 1 UT . TRUE and 


1 


NNEP INPUTS 

thermodynamic' 

INPUT 1 ,1 . 


START 


NNEP 

THERMODYNAMIC 

PROGRAM 


PERFORMANCE 

DATA 

ARRAY 



Figure 31. Overall Program Structure 
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IPLT = TRUE signals that the weight and dimensions calculations are to be 
performed, and the WATE-2 Data Set is required following the NNEP 
inputs. Figure 32 shows a typical card-stacking arrangement necessary to 
operate NNEP and WATE-2. 


2.2 Inout Description and Format 


The WATE-2 inputs are free-field format (NAMELIST) , and begin in Column 
There is no specified order to the inputs; however, for the following 
discussion they have been grouped into Plot-Print Indicators, Length 
Indicators, Mechanical Design Indicators, and Design Values. Figure 33 
shows a complete input set for a typical case. 


2.2.1 Plot-Print Indicators 


IWT 


IPLT 




= 0 - Do not do weight calculation 


1 - Turn on the thermodynamic parameter maximization of the 

WATE code. Do not do weight calculation 


? - Do weight calculation using maximum thermodynamic 

parameters 


3 - Do weight calculation but do not write maximum conditions 

for the components 


4 - Do weight calculation with airflow scaling 


= T - Gas path lavout 


F - No gas path layout 


placed at end of NNEP inputs to signal beginning of 
WATE-’ inputs 
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f 



b 



Figure 32. Deck Stacking 
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INPUTS ADDED 


t 




MOLmF 1 MOW BEING U3 ED 

Lm 

I P LT =T t 
1 3 11 =F , 

131w=F, 

I JU1CJ=4. • 

I L tMG ( 1 ! =2 ,3,3, 6, 7, 8, V, 10, 11* 

**MfcC( 1 ,2)=*FAN *, 1,1, 4,3*1, 

I wMtC( l ,31= 'SPIT*, 6*0, 

1 w ME C ( 1,4) = »iaJCT',C,5*0, 

I * Mt C < 1 , 3 1 = • lP *, *,l,2,4*c, 
i m MeC 5 l ,6) Pb UR * * 1* 5*0* 

IWMtwl l,7) = »r*>I •* 0*5, -5,3*0, 
l « Mt C ( 1 » d ) = * LP T ' » 1, 2, 7* 3* ( » 
iaKECU,*)*MIX ' * o*0» 

I*MEC< a » av ) = *JUCT • ,2,4*0, 

I « Me 2 ( 1, li I = *N 0 Z *,1,10,4* „, 

1 *MSC ( 1 « 12 »*»SrtAF» ,1,8 ,3*0 ,2, 

1 a Mt 1, a 3 > - *Sn4r • ,2,7 ,3*0, 5, 

0 - ^ V a L l a ,2 ) — *524,1 *7 , , 45 ,1 .5, 3*5, 4*5 *• 45 ,C.,0.,l.,t»,2*,l., 
Jt SVAL( 1 , 3 ) = 1>*0. , 

Je 3VAe( a ,**) = .45,2*0.,i 1. ,1 1*0. • 

Jt. «>VAl1 1 ,5 I = .4 3 , 1. 35 , . 70 , 1 .2 ,2 . , 1.3, .3,0. , 0 . , 1 .,0. ,2 . , 1 . , 

J c 3V At. I 1 ,4 ) = it)0.,.015, 

Jt 37 AL i a , 1 1 * .3 , .25 , A .3 ,1 .5 ,1.5 , • 53,1 30 LoC • , 3 . , 1. ,6 *0 . , 

J — 3V At ( i ,2 ) *• 3 5 , • 2 43 , 1 .5 , 2 . , 5* ,.b,15CCC0.,^.,l . , b * C . , 

Dz 3 Val < 1,4 l=l>*0., 

J u 3V A L ( A , 1 J ) - .1,2. , 

Jt SVAl < 1 , i 1 ) =1 • , 14*0. , 

U ~ 3V At_( a ,121 -5000o . , .5 , • 65 ,12*0. , 
jw oVAL 1 i ,1 J ) = ju0jj ., .3 ,1 3* C. , 

L t. Mj 




Figure 33. W ATE-2 Input 
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I SI I 


T - SI units Input 


I SI 0 


IOUTCD 


I OUT CD 



F - English units Input 
T - SI units output 
F - English units 

0 - Short form-engine weight, length, and maximum radius 

1 - Long form-component weiqhts and dimensions and short form 

? - Debug option and long and short form 
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2.?*2 Length Indicators 


z^i:r ^ «»» owt. M ^rib*. ,„ «•» 

ILENG in the”order 'that^the* 1 * " NEP CmB °™ nt ninber '* Rifled in 

-.a, r'V" ,n ,en9th to ach( - '* 

furthest downstream nozzle. Figure 34 ho C “"7 SS ° r lnd f " ld wi "’ «* 
'«* inputs for that englni. r L « “ * lVBiCi, - ** 

°«" ,e or shafts (13, and (14, heca„t These M ’’ 

contrihute to the total en qi „e length COmB °" enU d ° " ot 



ENGINE LAYOUT 



ILENG (I) • 1 .7.3,6, 7.8.9. 10. 1 1 . 1 ?. 


FLOW PATH ANO COMPONENT NUMBERS 


Figure 34. Length Input 












?.?.3 Mechanical Design Indicators 


e mechanical design indicators lUMEC) must be specified for each 
component of the NNEP simulation, ,1th the exception of the H«P 
Controls .Inlet. and Mater Injection or any „ t „er CMponMt „ ot 
represented In WATE-?. 


A number of shaft c<-ponents may he required to simulate an engine In 

’ as shown ,n Figure 7S. MATE-? will determine the weloht only for 

conneain, shafts of major components, such as the typical HP or lp 

shaft 1„ the example of Flgur, 31, only shaft is and shaft 17 would be 

specified. The smaller component number must always be used on the inner 

s a , With increasing component numbers as concentric shafts are added 

around the Inner shaft. 

IUMEC is a two-dimensional Integer array that contains ,11 of the 
mechanical design indicators. It is of the form 1MHEC (N, HI. where M is 
the component number used in NNEP, and N Is the variable number as 
OC , nod below for each component. Each variable in the IMHEC array for 

format! >mPOn ' ,nt ” S> ' 0 "" F ' qur? ,3 Tree-field hAHEUST 



SH 


INNER SHAFT 



IWMEC (1,151 « - 

IWMEC (1,17) ■ • 


Figure 35. Shaft Input 
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2 . 2.3.1 Compressors 


IWMEC Array 
Location 


Description 


Tvpe of compressor being weighed. 

'FAN' - Typical fan 

' FO ' - Outer portion of non-rotating splitter fan 

'FI' - Inner portion of non-rotating splitter fan 

'RSFO' - Outer portion of rotating splitter fan 

'RSFI ' - Inner portion of rotating splitter fan 

'LPC' - Low pressure compressor 

'HPC' - High pressure compressor 

This indicates if the fan or compressor has stators or if 
the compressor is a centrifugal compressor. 

1 - Stator weight is calculated 

0 - Stator weight is not calculated 

2 - Centrifugal compressor 

This is the indicator for 'front' frames in compressors. 
This input mav be: 


- No frame 

- Sinole hearing frame for turbofans and 
turbojets with Power Takeoff (PTO) 

- Single bearing frame with PTO 

- Two bearing frame, such as the frame in 
front of the HPC in the ,1T8D or JTPD which 
extends outward to the fan outer case and 
holds two hearings with PTO 


V ■ 


Location 


4 




f 


7 


Description 


This is the indicator for the 'rear' frame in a compressor 

0 - No frame 

1 - Single bearing frame for turbofans and 

turbojet without Power Takeoff (PTO) 

2 - Single bearing frame with POT 

4 - Two bearing frame, such as the frame in 

front of the HPC in the JT8D or JT9D which 
extends outward to the fan outer case and 
holds two bearings with PTO 

This is the component number connecting to this component 
for split flow compressors only. If this is the Fan 
Outer, the Fan Inner must be specified. If this is the 
Rotatinq Splitter Outer, the inner splitter must be 

specified, and vice versa. 

Gear box indicator - 0 - No gear or component number of 

shaf- 

Number of stages 
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Compressor number from which the RPM is determined 

Component number from which the outer radius limit for the 
turbine is determined. If the component number Is 
positive, the outlet dimension is used. If negative, the 
inlet dimension is used. If 0 , it will use the outlet of 
the feeding component. 

Number of staqes 

Indicator for axial or radial turbine 



k 


2. 2.3.3 


Burners 


Location 

1 


2 

1 


Description 

This is the type of burner being weighed. The input is 
the burner name in four spaces. 

' PBUR ' - Primary burner (airframe will be included) 

‘DBUR 1 - Duct burner (a mean radius is specified! 

•AUG' - Augmentor (no inner wall) 

This is the indicator for frame weight, normally only for 
primary burners. This frame includes a bearing. 

0 - No frame 

1 - Frame 

Ducts 


Indicator as to type of duct 

1 - Dummy - i.e., no weight or length 

? - Length input 

? - Length derived as in a duct connecting a 

splitter and a mixer 

4 - Cross over duct for centrifuqal compressors 

5 - Diffuser for centrifugal compressors 
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2. 2. 3. 5 Shafts 



? 


» 


I 


» 


1 


t 


» 



Location 


Description 

1 

1 SHAF ' 

- Standard shaft 


2 Shaft number from Inner to outer, i.e., 1, 2, 3, 4, or 5 

3-6 Turbine numbers connected to this shaft. The last entry 

Is the furthest downstream turbine. This Is used for 
power summation. 

7 First upstream compressor connected to the shaft 

2.2. 3.6 Mixers 

1 Type of mixer 


'MIX' 

- The coannul ar emergence of 
without mechanical mixer 

two streams 

' FMI X 1 

- Forced mixer, mechanical, 
lobed mixer 

i.e., Daisy 

Indicator for primary Input node 


0 

- Primary Is inner 


1 

- Primary Is outer 
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2. 2.3. 7 Nozzles 


Location 


Description 


'NOZ' 


- Input 


Nozzle type 


- Convergent 

- C-D variable area 


Component number from which the nozzle inlet diameter can 
be determined. If this diameter is taken from the inlet 
of the component, the (-) component number must be 
entered. If (+), the exit node will be used. If the 

previous component determines the diameter, this location 
may be zero. 


Thrust reverser type 


0 

1 

2 


- None 

- Fan 

- Primary 


The calculated component weight can be adjusted by an imput scaler, 
DESVAL (IS, M), Which is a factor applied to the calculated weight. A 
zero value, however, denotes that no scaling is used. If it is desired 

to zero-out the weight of a component, the scaler can be set to a trivial 
quantity such as .0001. 



2. 2. 3. 8 Splitters 


Location 

1 

2 

2 . 2 . 3 . « 

1 

2 

3 

4 


Description 

•SPLT' - Input 

1 - If inner stream is not primary 

Annulus Inverting Valve 
Input "VAIV' 

Location of Valve 

1 - Inner 

2 - Outer 

Component Number of Opposite Duct 
0 if Fixed, 1 if Movable 
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2.?. 3. 10 Heat Exchangers 


I 


* 

) 

> 


i 


I 


Location 

1 Input 'HTEX ' 


Description 


? Type 

1 - Fixed tube 

2 - Rotary 

3 How Direction 


1 

? 


?.’.4 Design Values 


- Parallel flow 

- Counter flow 


This section contains the mechanical and aerodynamic design data 
necessary to determine the weiaht and dimensions of each component. A 
summary of this array is shown in Table 4. If desired, the default 
values. Table F, can be used for any component by not specifying the 
inputs for that component. The data required is in the floating-point 
two-dimensional array DESVAL IN, Ml, where M is the component number from 

NNEP and N is as defined below. A typical range of values is shown in 
Table fi. 


Design limits are built into the proqram, as shown in Table 7, and cannot 
be altered bv inputs. If these limits are exceeded, the calculation 
continues and a warning is orinted out.. 


f 


8 ? 


■^.POSITION 


Table 4 PESVAL/DEFAUL Array 


-^4 


COMP 

TURB 


BURN 

DUCTS 

TRAN/ 

SHAFTS 

MIXERS 

AIV 

HEATEX 

NOZ 

SPLT 


f> = 


1 

MNI 

MNI 


VR 

MACH 

STRESS 

L/H 

L/H 

#TUBE 

L/D 

MNI 


PRM 
TLP * 


TR 

L/H 

RHO 

NO. PASS 
NO. PASS 
MNIP 

H/T 


H/T 

SOLID 


I DIA MEAN 
| DIA MEAN 

H/T 

MNI 

MNIS 


4 

5 

6 

SOLID 

ARI 

ARO 

ARI 

ARO 

MNO 

refloc 



refloc 



MNO 

RH 

WTIC 

BPR 

. 



MNO 

REFSTR 
.2% YIFLD 
STRESS FOR 
DISK 


WTOC 


2 9Ja 


h/N 


STAGES 




13 


RPMSC 


14 


TMET 


15 

WEIGHT 

SCALER 
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Table 5 DEFAUL Amy 
























TABLE 7 DESLIM ARRAY DEFAUL TYPE AND VALUES 



POSITION 



BLADE PULL STRESS CAN NOT EXCEED: 


FAN AND COMPRESSOR: 
!)P TURBINE: 

LP TURBINE: 


80000 PSI 
50000 PSI 
60000 PSI 


HUB/TIP FOR ALL COMPRESSORS CAN NOT EXCEED: 


HUB/TIP CAN NOT BE LESS THAN: 


FAN AND COMPRESSOR: 
TURBINE: 


0.50 


TURBINE STAGE LOADING INPUT CAN NOT BE LESS THAN : 
TURBINE: 0.28 


FIRST STAGE ALLOWABLE PRESSURE RATIO CAN NOT EXCEED: 


COMPRESSOR: 


LAST STAGE EXIT MACH NUMBER CAN NOT EXCEED : 
FAN AND COMPRESSOR: 0.6 


BLADE HEIGHT CAN NOT BE LESS THAN: 


COMPRESSOR: 


0.4 INCH 


NOT USED 
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2. 2. 4.1 Compressor 


Array 


Location 

Description 



1 

Compressor face inlet Mach number 



2 

Maximum first stage pressure ratio 



3 

Compressor face hub-tip ratio, R h /R t 



4 

Blade solidity, rotio of blade cord to 1 

blade spacing 

5 

Blade aspect ratio at first stage 



6 

Blade aspect ratio at last stage 



7 

Compressor exit Mach number 



8 

Maximum compressor inlet temperature. 

ZERO 

if design 


point temperature is to be used for 
°R, °K. 

material 

selection 

Q 

Maximum compressor outlet temperature. 

ZERO 

if desired 


point temperature Is to be usod for 
°R, °K. 

materai 1 

selection 

10 

Maximum speed ratio - RPM /RPM J . 

max design 



11 

Blade material density. ZERO if WATE-? is 

to select 


material. Ib/in , Kq/cc 



1? 

Compressor design type 
1. Constant hub radius design 
7. Constant mean radius design 
3. Constant tip radius design 
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Location Description 

13 RPM, scaler, normal input is 1. - use to match known RPM 
of engine 

14 Temperature at which a change of material is required. If 

ZERO 1160°R will be used, °R, °K. 

15 Compressor weight scaler, input ZERO if no scaling is 

desired 

16 Stator blade taper ratio. ZERO input sets 1.8 for fans; 

1.2 for compressors 

17 Blade volume ratio. ZERO input sets 0.055 for fans; 0.1? 

f«>r compressors 

Centrifuqal Compressors 

Location Description 

1 Compressor inlet face Mach number 

Maximum first staqe pressure ratio 
'< Compressor hub tip ratio 

a Hf’M ratio 

6 Compressor exit Mach number 

6 Clear ratio of the power shaft 

1 Horse power of power shaft 


ft- 17 Not used 
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1 


w 


Turbines 


Location 

1 

2 


3 

4 


8 

7 

•3 


n 




Description 

Turbine face inlet Mach number 
Turbine loading parameter 


U T ? /*rW" sti!| es 

Blade solidity, blade cord/blade spacing 


Blade aspect ratio of first stage 


Blade aspect ratio o f last stage 


Turbine exit Mach number 


OiSk refw-CTCf Stress . .■>* ^, {n \ 

turbine design tvpe 

1. Constant tip '•adius design 
7 . Constant mean radius design 
Constant bub radius design 

Maximum speed ratio - RPM /RPM 

max design 

Turbine control radius inches/cm - blank if transferred 
f rom a component 


Dens t tv of material in turbine blades - lb/in^/Kg/cc 

Blade volume factor. 7ER0 input sets 0.1SS for high and 
intermediate turbines; 0.10S for low turbines 

Not used 
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Location 


Description 






» 


I 


» 




♦ 


i5 Turbine weight scaler, input ZERO. If no scaling is 

desired 

lb Turbine blade taper ratio. ZERO input sets 1.0 for all 

turbines 

17 Stator blade volume factor. ZERO input sets 0.155 for 

high and intermediate turbines; 0.195 for low turbines 

Centrifugal Turbines 

1 Turbine face inlet Mach number 

?-5 Not used 

fi Turbine exit Mach number 

7-17 Not used 

?.?.4. 3 Burners 


Location 

1 

2 

3 

4 

5 


Description 

Burner through-flow velocity, ft/sec, m/sec. 

Burner airflow residency time, sec. 

Berner mean diameter, in. or cm. If zero, diameter is 
calculated to match connecting component 

Component number for calculating mean burner diameter. 
Enter zero if diameter is specified 

Number of cans for can burners 
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Location 


Description 


6-14 Not used 

15 Burner weight scaler, enter ZERO. If no scaling is desired 

16-17 Not used 

2. 2.4.4 Ducts 

1 Duct Mach number 

2 Length to height ratio of duct, required if mode 2 is used 
in IWMEC 

3 Duct mean diameter, in. or cm. If 0., duct diameter is 
calculated based on node specified below 

4 Node number to calculate mean diameter. Enter 0, if mean 
diameter is specified. Enter -1, if connecting component 
is to be used 

5-14 Not used 

15 Weight scaler, ZERO. If no scaling Is desired 

16-17 Not used 
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2.2.4. 5 Shafts 


Location Description 

A A 

1 Shaft allowable stress. lb/in\ Newton's/cm' 

3 

2 Shaft material density. Ib/in , Kg/cc 

3 Diameter ratio of shaft °^ nner / D 0 vjter 

4_7 Component numbers for total spool inertia 

8-14 Not used 

15 Shaft weight scaler. ZERO if no scaling desired 

16-17 Not used 

2. 2. 4. 6 Mixers 

1 Effective lenoth to diameter ratio of mechanical mixer, 

L/^2A/*, where L is the mixer length inlet to exit, A is 
the total flow area. Enter 0. if not a mechanical 
(forced) mixer 

? Number of passages 'or lobes) in mixer of either hot or 

cold stream. 

3-i 4 Not used 

15 Weight scaler. Enter ?ER0. !f no seal inn is used 

16-17 Not used 
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?.2.4.7 Nozzles 


Location 

1 

? 

3-14 
15 
1 6-17 
*.?.4.R 
1 

? 

7-14 

15 

16-17 


Description 

Length to diameter ratio of nozzle 

Bypass ratio for mixed flow nozzle for T/R weight 
Not used 

Weight scaler. ZERO. If no scaling desired 
Not used 
Spl i tte^s 

Only input if first calculated component in flow path Mach 
number in. 

H/T ratio in. 

B1 ank 

Weiqht scaler 
Not used 
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2.2. 4.9 Annulus Inverting Valve 


Location Description 

1 Specific length - L = V^A/n 

2 Number of passages. 

3 Mach number of inner passage. 

4 Mach number of outer passage. 

5 Hub radius in inches/cm or - component number from which 
hub radius is taken or blank if feeding component 
determines the hub radius. 

6 Inner cylinder weight - lb/ft?, Kg/M?. 

7 Outer cylinder weight - lb/ft?. Kg/m?. 

8 Wall weight - lb/ft^, Kg/M?. 


Blank. 


Weight scaler. 


16-17 


Not used. 


2.2.4.10 Heat Exchangers 


Number of tubes if "Fixed" type. 


Mach number in primary stream. 


W 


3 


Mach number in secondary stream. 


Engine Bypass ratio if "Rotary" type. 

Not used. 

2 . 2 . 4. 11 Miscellaneous 

"ACCS" is a one-dimensional namelist array that contains the value of the 
accessory weight scaler. Default value is 0.1. 

"DESLIM" is a one-dimensional namelist array that contains the mechanical 
design limit values for the components. It can have 15 values. First 13 
values are defaulted. Range of values is shown in Table 4. 

"I SCALE" is a one-dimensional namelist integer array which controls the 
engine scaling logic of the program. 

Output indicator 

Debug option and long and short form for every scaled 
engine point. 

Debug option and long and short form for unsealed engine. 
Long form for each of the scaled engines. 

Number of scaling points default is three. 

Not used. 


one-dimensional namelist array that contains values of 

scalmq factors. It can have six values. First three values are 
defaulted to 1., .8, 1.2. 


iscale(^ 

1 

? 

I SCALE {? ) 

I SCALE f 3 ) 
"SCALE" is a 


4 

5-17 
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"ACCARM" is a namelist array that contains the value of the centroid 
distance for the accessories component in the CG calculations. If no 

value is input, accessories are not included in center of gravity 

calculations. 

"DISKWI" is a namelist array that is used as an indicator for the new 
disk weight method. 

0 Do disk weight calculations using the old method. 

1 Do disk weight calculations using the new method. 

2.3 Program Output 

The output from WATE-2 may be selected in any of three output formats. 
Either English or SI units can be selected. Examples of the output are 
shown for the short output in Figure 36, the long form. Figure 37, and 
the debug output. Figure 38. This output shows the mechanical design and 

weight breakdown within the individual component. The units in the 

output section are shown in Table 8 for English and SI units. The type 
of units used are noted in the units section of the output. 

A flow path layout is also available for conventional type engines. A 
typical layout is shown in Figure 30. The layout is scaled such that it 
will fit on one page of the output. 
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lLjkAL L-Ant tNbl.Mt wclGHT* 2 91 i,. AcLtSSOKltS* 2 62 .0 ^ 
tsr Hated TjT AL LENGTH* 2J6. ESTIMATED MAXIMUM KAD1US = 


Figure 36. Short Output 


WtlbHl INPUT DATA IN ENGL JN1TS 


*£ l 

uh r u 

UT PUT DATA IN 

ENGL 

UNITS 






wT 

COMP 

A CCU 

UPSTREAM RADIUS 

DOWNSTR 

NJ 

cbT 

LtM 

LEN 

KI 

RU 

R1 

RO 

R I 

RO 

1 

c. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


Iu2v. 

29. 

29 , 

9. 

2 u) • 

c. 

0. 

13 . 

18. 


0 • 

c . 

29. 

C. 

0. 

0. 

0. 

IS. 

16. 

*♦ 

29. 

62. 

91. 

16. 

Id. 

0. 

0. 

16. 

18. 

5 

Cl 6 * 

26. 

59. 

9. 

13. 

(. 

0 . 

11 . 

11. 

6 

250. 

lb. 

72 . 

9. 

13. 

0. 

G. 

9. 

13. 

1 

126. 

6 . 

78. 

10. 

* A. 

0. 

0. 

10. 

13. 

3 

9o9. 

13. 

91 . 

9. 

12. 

0. 

0. 

9. 

14. 

9 

0. 

3 . 

91 . 

9. 

lo. 

16. 

21. 

A. 

21. 

10 

<♦6. 

63 . 

199 - 

0. 

29. 

0. 

0. 

0. 

29. 

11 

2v5. 

5b. 

2C6. 

C. 

2 9. 

C. 

c • 

G. 

27. 

12 

*♦0. 

0. 

0. 

9. 

20. 

10. 

11. 

0. 

0. 

1 j 

U. 

0 . 

0. 

9. 

13. 

0. 

0. 

0 . 

0. 

ILjIAI 

l An t 

tHGINE WEIGHT* 

2 911,. AC Lb SSOKIES- 

262 .0 2 


RADIUS 


R 1 


RU 


NSTAGE 


0. 

0. 

16. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


0. 

0. 

18. 

0. 

0 . 

c. 

0. 

0. 

0. 

0. 

0. 

0. 

0 . 


ESTIMATED TuTAl LENGTH* 2 36. ESTIMATED MAXIMUM RADIUS* 29 . 


Figure 37. Long Output 
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OOOOOfsJN>OOOOCOO 


» 


» 


I 


I 


t 


t 


T 


************* 

* * 

* t-AN 2 * 

* * 

*444**4*44**1 

oucr 

M NQ vtL T TOT P TOT P STAT AREA GAM 

0.524 570. 514. 1405. 1574. 6.9517 1.400 5 

U TIP STRESS OEN W/AREA TR H/T 

1258.9 26757.6 0.168 2.339 1.800 C.450 

COMPRESSOR 2 MECHANICAL DESIGN 

UJAD1NG N STG OIAK D TIP C RPM C RPM 

0.574 3.00 39.98 1258.9 7216.9 7216.9 

FRAME WT = 95.67 

ST AGE 1 

WD Kb WS WN WC CL RHOB RHOD AR 

65. 59. 59. 0. 26. 7.4 0.168 0.168 3.50 

PK DEL H MACH AREA K HUB R IIP NB U TIP STR WEIGHT 

1.4789 lo .7 0.524 6.95 2 8 . 99 1 9. 44 5 9 1 2 5 8.9 26 75P. 209. 

STAGE 2 

WO W6 WS MN WC CL RhOB RHOD AR 

91. 3h. 34. 51 . 21. 6.2 0.168 0.168 3. CO 

PK otL H MACH ARfcA K HUB R TIP N8 U TIP STR WEIGHT 

1.4155 16.7 0.499 5. 180 11.02 lb. 96 67 1193.2 20191. 231. 

STAGE 3 

WD wB WS WN WC CL RhOb RHOD AR 

V 7. 23. 23. 46 . 19. 5.7 0*168 0.168 2.50 

PR UEL H MACH AREA R HU8 R TIP NB L TIP STR WEIGHT 

1.3o71 16.7 0.476 4.017 12. 17 18.23 70 1148.1 15768. 208. 

FRAME WT = 2 r . 5.1 5 


t 


t 


N STG wElbHT LENGTH 
3 1C23 .68 28.80 


DOLT 

M NO VtL T TOT P TOT P STAT AREA GAM 
0.*t50 562. 727. 5447. 4743. 3.2206 1.3951 


PR AD EF 
2 .8600 U.B700 

Hi HC 

123.9 5 174.07 


PU TO HP 

5447.2 726.9 16910. 

W1 CWI 

238.50 265.00 


******************* TOTAL COMP WEIGHT IS 1028.680 

Figure 38. Debug Output 
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UJCl 

M NO VEL T IOT 
0 . 9* 0 5b2 • i £ 1 • 


P TOT P STAl AREA GAM 
5W. *J5 3. 1.81*6 1.3951 


U TIP 5TRESS 
liibS. 1 23 3-'' 1.5 


OEN W/ARcA 
0.168 C.687 


TR 

1.20C 


H/T 
C. 700 


CUMPRtbSJK 5 MECHANICAL Ufc 51 ON 


L CAUl Nb 
o • t) 5 1 


N S fo 

1 U# ou 


OIAM J TIP C R*>M G RPM 
25.58 1085 .6 115 15 .5 9727. > 


FRAME WT = 118.22 

1 


STAoE 1 CL RHUB RH JO AR 

* d .f 3 ? l£ 9.5 0.168 0.168 2.00 

15. 15. 36 ‘ TI p STR WEIGHT TIN 

•* X a * 1 9*. 727. 


S, ‘" E ^ WS WN WC CL Rh06 KHOO AR 

-°- *• *• 2 *‘ «• °; l ? t ip 0, Nt 8 U iir STR WEIGHT TIN 

25z h U-!^ 6C 125C.7 18516. 75. 600. 


s l* ;: 4 :! 2 ;: f. c m 0 ^s /.s, 

KK JLL H MACfl AREA A HU8 K . I 1 L *7^1225.3 15058. 59. 873. 


WEIGHT TIN 


rw w l. w ^ \ m 

l..:9to 17.o C.h 2C I.l7i 9 . ? 5 12.19 


-o' WS -N -C CL Hh3» K»3D ** 

) i 4. ? >b \j . G • 1 ob 1 • 83 

lj 2 *1 H «a K HUS K UP NB U TIP STR -EIGHT TIN 

iMu, L7.6 e.,Cu 0.970 9.9b 12.00 01 1206.2 12*77. *9. 9*6. 


MiT -o' -S -N -C CL KH0» KHOt *« 

“• 3 ' J ; K hGs i UP mS U UP STR -EIGHT TIN 

"Mo Mil. .»>.» . 0 = 27 . *>. 


S ’M -8° WS -N -c CL KHOB RHOO AR 

‘pk DEL H ^PAlH° *AR t A* K HOB MtIpMS U TIP SIR ,<E ’£ HT 

17.6 0.37a c.701 10.28 11.75 105 1180.1 9C18. 36. 108 . 

Figure 38 . Cont. 
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WL> 

W8 

WS WN WC CL KHUB KhOD 

AR 

V. 

4. • 

2. It. 4. 1.8 0.168 0.168 

1.67 

Hk 

1 1 3 » 

utL H 

MACH ARtA K HUB R TIP Mb U 

TIP 


1...1/1 i/.o 3.360 0. 60 3 10.39 11. 65 116 1171.0 7829. 

S i Abt d 

Wu wt WN WC CL RNUB RhOO AR 

ic * ->• J » 13. 1.6 u.2Bo 0.286 1.61 

RK uu H Ma«-H AKtA K HUB k HP NB U TIP SIR 
1 . 17.0 0.3t5 0.534 10 .A 7 11.68 126 1163.6 11712. 


ST AG t 9 

WD wb ws WH WC CL 

lb. j> . 5. 12 • 5. 1.5 

«-K utl H MACH ARtA R HUB 
1 17.6 *,.330 0.0/6 10.6 


KHUB RHOO AR 
0.286 0.286 1.66 
R IIP NB U TIP 
11.62 13 6 116 7.7 


STR 
10 AO 7. 


Si Abb 10 

wt > ** WC CL KHUB KrtOL AR 

1 IT • 3. 1 .A 0.2b6 0.286 1.60 

HK uul h MAwH ARtA R HUB K TIP N 8 U TIP STR 

l.lb^H 17.o 0.315 O.A2t 10.5V 11. 47 1**6 1152.8 °3A3. 


N STb WEIGHT LtNGTh 

10 616.43 25 . A 3 

dolt 

* V6L 1 r3T P WT •» STaT AKtA GAM 

O.ibb 1 A 36 . 5123 o. 402 ^i. 0.3874 1 . 353 V 

PK Au fc> PU T ( HP 

^.V 06 U O.b/ob 5 i 235 • 9 1435.6 3 JV 65 . 

Hi HU *1 GW 1 

l 7v . 0 7 3 52.2 3 l z>4 . 7 5 61.9 7 


total comp wtiGni is 


616 ,v77 


* * 
* 


♦ 


* PBJK o * 

* * 

** *♦♦***7 

O'ORNti NUMBtK o 

kin rout llngth 

6. / 5o 1 2.909 IB .000 
CAS WT LIN Wl NO l W.T 

2-*.2 40. v 17.8 


MACH 
0.055 
INC WT 
1 O . 4 


WSPLC 
4 • 59t> 
F-KAMfc 
15 1.3 


WT ITT 

26 0.2 


Figure 38. Cont. 
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♦ » ♦♦♦ * *♦+* **2 

l)U CT 
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0.5U0 UEo. 2o2i. 90ii2. 


p SI A I AH £A GAM 
393 21. 3. ->977 1.2 9o8 


o UP STRESS OtN 

1 luo. 0 9819.5 0. 2b 6 
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OEL H RPM TGRG 
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SI Aot 
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0.5 ->o 1 1 *»9 . ^.628* 1 29 36 . 
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H IN 
o vv. c a 
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H 00T 
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AREA 
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jj 9 6 V. 
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2027.7 
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TOTAL T0R8 wtlGHT IS 
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* lp r c * 

* * 


DJLT 

M NO VEL T TOT 
0.560 11*9. 2vi2U. 


P TOT P STAT 
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AktA CAM 
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l 
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Figure 38. Cont. 
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Figure 38. Corn. 
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WEiGHT INPUT UATA IN ENGL UNIT* 
WEIbhl OUTPUT UATA IN ENGL UNITS 
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Nu 
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0. 
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t»7a. 
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*16. 
0 • 
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Figure 38. Cont. 
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2.4 Sample Cases 
2.4.1 Large Engine 

A simple mixed flow augmented turbofan is used as an example for the 
WATE-? input and execution. Figure 40 shows a schematic and a block 

diagram of the engine. From this block diagram, the component numbers 

are determined. 


To construct the input deck the indicator section must first be set, 
Figure 41. In this example, the units in and out are English, so -SIT 
and I SIO are set false. Since the weight and gas path Uyout 'are 
desired, IWT and IPLT are set true. The debug option is turned on with 
TOUTED set equal to 2. The length inputs are then entered in ILENG. 
Since the duct (4) and shaft (12) and (13) do not contribute to the total 

length, they are not entered. Also, the components are entered as the 
flow would progress through the engine. 


The IWMEC values are now entered. Since no inlet weight calculations are 
done, the inlet is not entered. This is true with any component entered 
m the NNEP KONFIG section; it is not entered in IWMEC if no routine 
exists to weigh it. In the example, the IWMEC (1,2) card says a "fan" is 
oemg weiqhed. The weight will include stators, IWMEC (2,2)=1, a front 
frame, IWMEC (3,?)=t, and an intermediate frame, IWMEC (4 2)=4 The 
IWMEC (1,8) card says a "LPT" is being weighed. It has a turbine exit 
frame, IWMEC (2,8)-l, and it is connected to component 2, (IWMEC 
(3,8)='/. The nozzle has variable area capability, IWMEC (2,11)=2, and 
its diameter will be taken from the inlet to the augmentor, IWMEC 

(3, 11, =10. Since the augmentor has constant diameter, the node position 
tor taking the diameter is of no consequence. 


^e DESVAL inputs follow the IWMEC inputs. Component numbers used in 
OESVAL must agree with those used in IWMEC. Input o, DESVAL data will 
override the default values. For the example case, the fan design card 
DFOVAL (1,2) indicates that the compressor inlet Mach number is 0.^24 
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MODE 1 NOW BEING USED 
LW 

IPIT = T, 

I S 1I=F , 

ISIO=F. 

IOUTCD=2, 

ILENG( 1)=2,3,5,6,7,8,9,10,11, 

I W ME C ( 1,2)=’ FAN ‘,1,1, 4,3*0, 

1WMEC( 1,3)=»SPLT»,6*0, 

I WMEC ( 1 ,<*) = • DUCTS 3, 5*0, 

IwMEC( 1,5) = 'HPC • » 1, 2* 4*0, 

I W ME C ( 1*6) = 'PBUR •* I* 5*0* 

I WMEC ( 1,7) = ' HPT ',C, 5, -5,3*0, 

1 W MEC ( 1 » 8 ) = ' LP T », 1,2, 7,3*0, 

I WMEC ( 1 ,9 ) = • MI X * , 6*0, 

IWMEC( 1,10) = 'AUG ',6*0, 

IWMEC( 1 ,11)= *NOZ », 2, -10,4*0, 

I WMEC! 1 ,12)=*SHAF» ,1,8 ,3*0 ,2, 

I WMEC ( 1,13) = *SHAF *,2,7, 3*0 ,5, 

Ot SVAL (1,2 >=.5 24, 1.7,. 45,1 -5 ,3 .5 , 2 .5 * . 45 , 0., 0 . ,1 0. ,2 •• 1 . ,0 * 1 .1 , 

DESVAL( 1,31=14*0. ,1.1, 

Dc SVAL ( 1*4) = .4 5,2*0., 1 1. ,10*0. ,1. 1 , 

0 l. SVAL ( i , 5 ) = «*/5 * 1 • 33 , . 70 , 1 .2 ,2. ,1.3, *3,0.,0.,1.,0. ,2 .,1*, 0,1.1, 

Ot SVAL (1,6) = 100. ,.015,0. ,5. ,10*0 .,1.1, 

Ot SVAL (1, 7) = .5,. 28 ,1 .5,1 .5 ,1.5 ,. 55,1 50 000., 3. , 1. ,5*0., 1. 1, 

Ot SVAL (1,8) = .55,. 243, 1.5 ,2., 3. , .6 , 150000 . ,3. ,1 . , 5* 0. ,1 .1 , 

Ot SVAL (1,9)= 14*0., 1.1, 

DC SVAL (1,10) =250., .016,12*0. ,1 .1, 

DtSVALt 1,12) =50000., .3 ,.85 ,1 1*0. , 1 .1 , 

Dt SVAL ( 1,11) =1., 13*0., 1.1, 

OESVAL( 1,13) =50000. ,.3, 12*0. ,1 .1, 

LEND 


Figure 4 1. WA TE-2 Input Example 
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the maximum first stage pressure ratio is 1.7 and the inlet hub/tip ratio 
is 0.45. The compressor has a blade solidity of 1.5 with a first blade 
aspect ratio of 3.5. The last stage has an aspect ratio of ?.5 and an 
exit Mach number of 0.45. The inlet and exit temperatures calculated in 
NNEP will be used for disk material determination, DESVAL (8,2) and 
DESVAL (9,2) are 0.; the RPM ratio between maximum and design is 1.0. 
The blade material will also be chosen by the code because DESVAL (11,2) 
is 0. The design of the fan is a constant mean line since the mode, 
DESVAL (12,2) equals 2. Also, no speed scaling or weight seal i .g will be 
done since DESVAL (13,2) is 1.0 and DESVAL OS, 2) is 0. A material 
change temperature of 1160°R will be used since DESVAL (14,2) is 0. 

The HP turbine DESVAL (1,7) has an inlet Mach number of 0.5 and a turbine 
loading of 0.28. It has 1.5 solidity with inlet and exit blade aspect 
ratio equal to 1.5. The exit Mach number is 0.55. The disk material is 
a high strength super/alloy with a reference stress of 150,000 psi. A 
constant tip radius is used in the design, DESVAL (8,7)=3., and a speed 
ratio of 1. is specified for stress calculations. 

To end the inputs, a "&END" is entered. This will initiate execution of 
WATE-2. The output of WATE-2 is shown in Figure 42 for the example case. 
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************* 
* * 

* l-AN 2 * 

* % 


duct 

M NU vtL T TOT P TOT P STAT AREA GAM 
0.52^ 570. 510. 1905. 1579. 6.9517 1.9005 


U TIP STRESS 
1258.9 26757.6 


DEN W/AREA TR H/T 

0.168 2.339 1.800 C.950 


OUMPKESSoK 2 MECHANICAL DESIGN 


uGAOlNG N STG DIAM 0 TIP C RPM C RPM 
0.6 79 3.0 0 39 . 98 1258.9 72 16.9 72 16 .9 

FKaMc WT = 95.67 


SI Abe 

WD 
9 1. 

PK 

i • 4i 5:> 

$ 1 AbE 

wO 

vy. 

PK 

1.3o71 
FH AME 


Wb Wi> WN we CL RH 08 RHOD AR 

26. 7.9 0.168 0.168 3.50 

UfcL H MACH ARtA R HUB R TIP N8 U TIP STR WET fjjT ttm 

' *'>•’ <>-«2 8.9, ,9.99 59 ,255.9 »S,. ^ 5 ^. 

2 

wo W5 WN we CL Rh08 RH3D AR 

3t. 39. 51 * 2 1. 6.2 G. 1 08 0.168 3. CO 

UtL H MACrt ARtA k *Ub R TIP NB U TIP STR WET rwT ttm 

,6.7 0.999 5. ,8 0 , 1.02 ,6.99 tT liil* 2»I5l. "HR S3. 

3 

W 8 WS WN WC CL KHO 3 KriOO AR 

23. 9 t>. 19. 5.7 0. lob 0.168 2.50 

uEL H MACH ARtA R HUS R TIP MQ L TIP STR WETTMT TTM 

,6.7 0.97> 9.0,, ,2., 7 18.2, 70 llla.l ,57^8. "US 658. 

wi = 285.15 


N STG wfloiiT LENGTH 
5 x C 23 • 63 2 8.80 

30 CT 

" r ^ 3 VcL T T 'J T P TOT p stat area gam 

t21m 3H97. 9/93. 3.2206 1.3951 

PK AJ C F PU TO HP 

2.660 0 o.b7oO 5*97.2 72o.V 16V10. 

Hl MC W I CW I 

123.95 17*. 07 9.38.50 265. OC 

*****♦**♦*♦**,**„, TOTAL WUGfi| 1S l028%fc80 


Figure 42. W ATE-2 Output Example 


**•»**♦*•#*** 


* HPC 5 * 

* * 

************2 

UJCT 

M NO VhL T 10T P TOT P STAI AKfcA 
0»45u i>b2. 727. 54 47 . -4743. 1.8196 1.3951 

U TIP STRESS DEN W/AREA TR H/T 

l<_65.1 233->1.5 u.loB C.687 1.2CC C.700 

CUMPRLSSJR 5 MtCHANICAL Ufc 51 UN 

LDAOINu N SIo 01 AM 0 T IP C RPM L RPM 

lO.ao 25.58 1085 .6 11515 .5 9727.5 

ERAM5 WT = 118.22 

STAoE 1 

Wl) *3 WS WN WC CL KhiJH kHJD AR 
24. 14. 14. 3a * 10. 4.5 O.lao 0 .1 o 8 2. OG 

PR ufcL H MACH ARtA R HUB R TIP Nb u TIP S TR 
1.3oGj 17.t> O.*.3o 1.320 3.4a lc.7v bO 1285. 1 23331. 


SlAbE 2 

WG W<- WS WN WC CL KKUti KHUD AR 

-0. v. 4 . 24. 6. 3.6 0.168 0.16 8 1.94 

PK jlL H MACrt ARtA R HUB k TIP Nb U TIP STR 

1.3^41 17.6 J »h3j 1 • 44^. 1 2 • 45 of 125 C. 7 16 51 1 >. 

5 I AGE _> 

* l> ^ *\ WC CL khOa KHOD AR 

la. e>. a. 25 . 7. 3.0 O.loO 0.166 1.89 

pk JtL k MALH AREA >< HUB R TIP Nd U TIP STR 
i..V*.o 17. a C.h 2L 1.17i 9.75 12. 19 7( 1225. 3 15046. 

SlAUfc 4 

WU ws WN WC CL KhOa KHOD AR 

21 . 6 . 2 . a u.lad C.loa 1.83 

PK ut.L H MAcM AKlA K HU3 K UP Nb 0 TIP STR 

l.ctQi i7.6 C.hCj 0.970 9.98 12.00 61 1206.2 12477. 

T Ao L 1> 

WJ *B WS xN we CL KhLB KHJL AR 

1 U J ‘ ->• lo » 2.2 v).l6d u.168 1.78 

Pk OEi_ H MmLH AKtA k HUB K Tl P Nd L 1 I P STR 

1.2499 1 7.P w. jVO a. bib 10.15 11.8a 93 1191.5 10 52 7. 

S f Aoc o 

^ * rt -<3 WH WC CL khOB Km 0 U AR 

1C * 3 ‘ lo * 4. 2.0 0. lob 6.1 ta 1.72 

PK GEL H MA 6 H ARtA K HUb K Tl P Nb U TIP STR 

i.t.J*4 i/.o Omilj 6.701 10. 2d 11.74 1 u4 1180.1 9C1P. 


Figure 42. Cont. 
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STAGE 7 

WO WB WS WN WC CL RHOB RHOD AR 

V. 2. 2. 14. 4. 1.8 0.168 0.168 1.67 

m UtL H MACH AREA R HUB R TIP MB U TIP STR 

1.2173 17.8 J .3 6 0 0 . 60 8 10.39 11.6b 115 1171.0 7829. 

SI AGE 8 

WO WB WS WN WC CL RHOB RHOD AR 

it. 3. J. 13. 3. 1.6 0.286 0.286 1.61 

Pk UtL H MAoH ARtA R HUB R TIP MB U TIP STR 

17.6 0.34b 0.534 10.47 11.58 126 1163.6 11712. 

ST Abt 4 

WO wb WS WN WC CL RHOB RHOO AR 

lb. 3. 3. 12 . 3. 1 .5 0.286 0.286 1.56 

Bk 6tc H MACH ARtA R HU8 K TIP NB U TIP STR 

1.14*16 w. 6 3.330 0.475 10.53 11.52 136 1157.7 10407. 


» 

»•* 


I 


I 


* 


» 


SlAbE 10 


W 0 

Wb 

WS WN WC 

CL 

RHOB 

RHOG 

AR 

i 

di. » 

2 . 11. 3. 

1.4 

0.2 86 

0.286 

1. 50 

HK 

GtrL h 

MACH ARtA R 

HOB 

K TIP NB U 

TIP S TR 

1. lb^ 

17. 

8 0.315 0.426 

10 . 5 9 11. 

47 1 46 1152.8 9343. 

N STG 

WEIGHT LENGTH 





10 

616 

.43 2 5.43 





DOLT 







M NJ 

VEl 

T TOT P TOT 

P 

stat 

ARtA 

GAM 

0*300 


W36. 61236. 

432 

-3 1 . 

0.3874 

l .3539 

PR 

AU 

EE PU TO 

HP 




v .406 0 0 .6 7uO 51235.9 1435.6 33965. 

Hi HU W 1 CW1 

1 74.07 352.23 134.75 61.9 7 


*****#*****+******* TUlAL COMP WEIGHT IS 616.477 


♦ * 

* PBJK 6 * 

* * 

* ♦ *♦* * *♦****?. 

BUKNtl NUMBtr\ < 

KlN k 0UT 
6. /bo 12.904 

CAS WT LIN Wl 
2**.2 40. V 


LENGTH 
16 .0 00 
NUZ WT 

17.8 


MACH 

0.065 

INC WT 
1 6 . 4 


WSP EC 
4.59t> 
ERAPE 
151.3 


wt or 

250.2 


» 


Figure 42. Cant. 
113 




L. MWWW ! Wb .'IWJiUAIMm* WM'IPy. M > WW ; IUiWTIPWlH.''. 


WEIGHT TIN 
32. 1159. 


WEIGHT TIN 
39. 1229. 


WEIGHT TIN 
36. 1299. 


WEIGHT TIN 
33. 1367. 


I 


\ 

J 


I 


1 


P TOT P STAT AREA GAM 
96112. 3 93 4 . 7 • 0.3977 i.29o8 

u TIP STRESS OEN W/AREA TR H/T 

lluo.O 9819.5 0.266 C. 296 1.000 0.922 

T UK31 Nt 7 MtCHANICAL OtSlGN 
H/I N ST G LO AO IN G AREA 

0.92.' 2.00u 0.2b 0 0.39b 

Ul KT1P RHOB DEL H RPM TORQ 

UOo.O 11.0 10.1 17 *»• 6 11516.5 1 8tR 13. 


S 1 Aot 1 

DISK. bLAUt VANE HWO CASE AR 

o.*» 2.3 8.4 21.7 4.0 1.50 

Kk JtL H MACH AREA R HUB R TIP NB U T IP STR WEIGHT LENGTH 

1.3956 fc7.3 0.5 0 0 0 . 398 10.19 11.01 180 1106.0 9620. 42.84 2.02 

SI Aofc 2 

1SK bLAUt VANE HWD CASE AK 

10.3 o.h 23.7 35.5 6.6 1 . 60 

PK Ull H MACH AREA R HUB k TIP NB U TIP STR WEIGHT LENGTH 

L.LuOj b7.^ 3 . 5l5 0.6O6 10. 19 11.66 116 1160.9 16956. 63.19 3.2^ 

U STG LENGTH WEIGHT 
2 5.31 126.03 

DUCT 

rt NJ VLL T TOT P TOT P STAT AREA GAM 
0 • 6 jU lltv. 2J26. 12436. 10243. 1.2074 1 .312 7 

PK TK AO EF PO TO TL.l 

3 .? 0 o I 1.292 8 0 . 860 0 12 43 5.6 2 02 7. 7 2027.7 

rt IN H OUT AREA FLOW HP 

099.4.3 64 - 9 . T9 5.17 1^7.56 .0964. 

*•*•*****♦****♦♦♦♦* TOTAL TUR 8 WtlGrlT IS 126. 028 

♦ ♦ 

* LPf C * 

* * 

* » *•*»**♦* **2 

DJlI 

M NO VLL I TOT P TUT P SI AT AREA GAM 

0.660 1 . ^j 23. 129 38 . 102 *.5. 1 .20 /_ t .3 12 J 


Figure 42. Cent. 
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* * 

* HPT l * 

* * 
************ 2 

OUCT 

rt NO VLL T TOT 
0.600 llbo. 2621. 


r 




U TIP STRESS DEN W/AREA TR H/T 

UI.G 11708.5 0.286 0.77 7 1.000 0.765 


TUKBINt 6 MLCHAN1CAL 
H/T N SIG LOADING 

0.765 2.000 0. 24 3 

oT RTIP RHU8 

727. C 11.5 8.8 


DE SIGN 


AREA 


1.20 7 


DEL H 

RPM TORO 

86.9 

7216.9 147693. 


STALL 1 

UiSK bLADc VANk HWD CASt AR 

5.0 2 2 . bo .0 38 . V 9.9 2.00 

P,x ulL H MACH ARkA k HUB R TIP N8 U TIP STR WEIGHT 
l.^ooV -.3.4 0.550 1.207 8.83 11.34 60 727.0 11709. 142.18 


STALE 2 

Lt 1SR oL AOt VANE HWD CASfc AK 

o.v 2 ?.o 81.4 34 . G 9.3 3.00 

K«, UtL H MAlH ARkA R HU8 R TIP NB U T IP STR WEIGHT 
1 . 81^0 4->.4 V -.578 1.632 3 . d 3 l 2..>9 98 780.6 16019 . 159.21 

FRAME *T = 167.79 


n length weight 

2 1j.*»1 **69.18 

uglT 

4 NO VLL T TOT P TOT P ST AT AKcA GAM 
O.tCO 1 i.34 . 1722. 5394. 4436. 2.3313 1.3249 


1 R 

.2 2 _> o 

h IN 
5 i*. . 7 / 


TR 

1 .1779 
H OlT 

*♦ . 1 / .8 7 


AO ef 
j .dooo 
AREA 
16.8 0 


PO 

5393. 7 
FLOW 
137.5 6 


TO 

1721 . 5 
HP 

16912 . 


TO. 1 
172 1 . 5 


TOTAL TUR8 WEIGHT 


IS 


469.184 


**«***♦¥♦*♦** 
» ♦ 


» AUl 1’J * 

* * 


II O KN t S .'(LMbtK 10 

KIN *0JI LENoTH 

u.o 24. 120 46.000 

las hT LIN WT NUZ W1 

2 3.7 12L.L 2 96.1 


MACH 
0.1**3 
INC WT 

O.C 


WSPtC 
11 .899 
hi 0 I 

439.8 


Figure 42 Cont. 


LENGTH 

4.77 


LENGTH 
4. 17 
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************* 


************2 

NUZZLt 11 

Wt 1GHT = 568.95 L ENGT H= 


-♦8.239 Trt Wl~ 


************* 


DOC I 


Ov»CT , 4 

Kn= 16.78 RT- 17.69 LENG= 62.16 

Ak Ea= 1.4CI RH0=.l6d 

CA6 WT INC WT WTOT 

l>.5**Ob 13.8 556 29.396** 


S HA t 12 


SHAhT 


D1 LENG 
3.81 98 . 74 


ON WT 

0.65 40.03 


brtA H- 13 


irlAhT 13 

L»J Ul LENG 

**.35 3.9** 18. CC 


ON WT 

1.27 14.33 




* ACCb WT 


*♦***♦ 


AwLS WT- _>ul.4l4 


Figure 42. Cont. 
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*=” # UN «I W S*" *}°«» OCIXNSTREAH RAOIUS 

* u. 0. „ . . RI R0 " W WTAS* 


* 1132. 

3 0. 

*♦ 32. 

i> o7d. 
t> 2/3. 
7 139 . 

a 316. 

7 0. 

10 <*a<*. 

1 1 C»_ o • 

12 9 h. 

*3 16. 


0 . 0 . 

9 . 20 . 

0. 

^ - 18 . 

'• 13. 

9. 13. 


0. 0. 

0 . 0 . 


*8 . 139. 

*&• 187. 


'• C- C . 

lo* P* C * 

10 « il* 0 . o. 

9 * lc« 0. 0. 

J[* i 6, *&• 21 . 

?• 0. o. 

°" (. c . 
9» 20. io. n. 

V - 13. o. o. 


°* 0. 0. 

18. o. o. 

16. i6. ia # 

18. o. 0. 

11 * 0 . 0 . 

13 . 0 . 0 . 

13* 0. 0. 

1^» 0. o. 

21 • 0 . o. 

0 . o. 
22 . o. o. 

0. 0 . o. 

8* 0. o. 


Tt,IAL UArtt ENGINE NE 1 GHT= ,g. , ’ °* °‘ 

e stihated tut At UE NGm- r, ,tUSS0R,ES = »«•*« 

LENGTH- 1 . 7 . ESTIMATED M4XIMU(1 rjc>jusi ^ 


42. Cont. 
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?.4.2 Small Engine 

A typical turboshaft engine as illustrated in Figure 22 has been 
represented by the inputs shown in Figure 43. The general input format 
for small engines is the same as for large engines. Major differences 
are the IWMEC inputs that call for centrifugal stages instead of axial, 
and the IWMEC input that calls for the optional disc weight procedure for 
the axial-flow turbine. 

Component input data should be selected to be representative of this type 
of engine. For example, blade volume factors may be significantly larger 
because of the small -size blades. Blade geometry also may be 
considerably different than the default values. Samll engines generally 
have larger blade taper than large engines. 

Output data for the sample case is shown in Figure 44. Since the engine 
geometry and component arrangement does not generalize for small engines, 
it is not possible to automate a computer drawn flow path layout as was 
done for the large engines. Radial component dimensional data also is 
not readily defined, except for the major envelope dimensions which are 
shown in the output. The same outDut options also are available for the 
small engines. 


2.5 Program Diagnostics 

The WATE-2 program contains error printout to aid the user in trouble 
shooting an input deck. A listing of the error messages and their 
meanings are shown in Table 9. None of these errors will cause 
termination of the program. The component routine in which the error 
occurred will be terminated and the program will continue its calcula- 
tions. The components calculated after an error may or may not be in 
error . 
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j 

j 

j 


su 

I «T=2,IP.TrT, • - - 

ISII=F,ISI0=F,I0UTCD=2, 

I L EM S < 1 ) = 2 * 3 * ft * 

3ISIWl = l., — - 

ft CCS = .1 7 

I U ME C < 1 , 2 ) = * HPC »,2,0,0,0,1,Q, 

IWMEC<l«3) = *P3UR*»L'«5*r* 

**5*?*2*3*i/*G* 

I y ME C < i ,5 ) = *M02 • * 1,5*0, 

lWMEC<l,b> = »SHAF»*l,A*3*'!*2, 

3ESVAL<l»2)=»7 £ *»b.*.ftF,.72Sb»,3,13.»143u.*l’3*0»» 

3tSVAL<l,3) = 3'J.*. 037*12. 9*2, *13*0*, 

5ESVAL(1,A) = «2,»A2*1*2,*9,3«1, •3**50 0 GO., 3«,1«, £>*(!. *3., 
DESVAL<1*5)=1.»1A*C., 

3 l S V A L ( i,t»)~tj"ji,v«*.3*«H3^2**A»*l„*i*| 
i E M3 


1.5 



Figun 43. Smail Engirt t inputs 
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* 

* CCNC 2 

* 


* 

* 

* 




*UX CDNiDiT IONS OCCUR AT 

**************************************** 


ALT 

MN 

VALUE 


p r o T : ■ . 

j • o v r 

1 A • 7 

Ls /SU1 N 

T T 0 T C . 

.( • 5 l b 

518.7 

deg r 

Z d I N 

r, • Tot 

*2.0 

Ls/SEC 


* * ****** + *****************★* -^ ****** ***** 


DJCT 

vj m 3 tfEi_ T TOT P TGT P STAT AREA GAM 
« ? A C 7 8 A * 519. 2119. 1A71. .259A 1.AM5 






1 

I 


I 

l 


S, 


B 


ST AGE 


ST PR 
3.5 30 

R ? M 
552 . b 

RT IP 
5 • 7 

INERT 1 a 
1 ft li • 9 

LENGTH 

3.5 

WEIGHT 
13. 8 

STAGE 
ST PR 
3.3 i j 

cl 

KPM 
5 D rT * O 

RT IP 
0 • 5 

INLRT I A 
91.2 

LF NGTH 
3.5 

UE Ii>HT 
1 3 . ■: 

* * * C j M P 

H3USI Si G 

* * * 





4 E I 3 H T = 2 1.7 

♦•♦CROSSOVER duct **» 
uiNjT J: . 3 ‘ 2 bE!SHT= 15.3 


**OI r F JSER*** 

L E R G T H = » 1 ■ a *EIGHT= 5 • S 

« * * u: M 3 ) X « * * 

SRP = 1A3«. GEAR K= 15. j SHAFT RPM = 2.2A. 

« E I GHT = 2.5.9 

^ STj b £ 1 5 ri 1 L £ N G T H CiNGRm I N _ R T I A 

2 id.ft 7.3A 3.7 192.1 

3 JCT 

vi |. w V t. T TCI P TGT P STAT AREA GAM 
,5]j 4‘i. lbi3fl. lby2D. •*6*0 1**7M 


HR AD IF 

3 ' 1 ' - 

HI HJ 

12 3.95 2 5 ! • • > 


PO 

1 3 9 <: D • 8 

. 1 

12 .‘ ' 


TO HP 

1.38.2 2 1 A ■•' . 

c:. 1 

I2.cn 


Figure 44. Srm/i Engine Output 
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* " * 

* PBUR 3 * 

★ * 

************2 


MAX CONDITIONS OCCUR AT 

********** ************ * * ********* 

* * * * * * * 

PT OT 

ALT 

0. 

MN 

P • £ L L 

VALUE 

117.6 

L3 /SuIN 

TTOT 

J » 


1936.2' 

DEG R 

CwIN 

a • 

0 • fl 0 v> 

2.1 

L3/SEC 


**************************************** 


Borneo number 3 





RIN ROJT 

LEN jTH 

MUCH 

WSPEC 


3.421 6.45b 

15.32 

.tu 

1.625 


CAS WT LIN WT 

NOZ WT 

INC WT 

FRAME 

WTOT 

2.5 16. 4 

8.6 

1.1 

0.0 

28.6 


************* 





* * 





* HPT 4 * 





* * 





************ 2 





MAX CONDITIONS OCCUR AT 




***************** 

********************* * * 


A_T 

MN 

VALUE 



PTOT 0. 

0 . C U c 

112.2 

L3/S0IN 


TTOT 


2100. C 

3"G R 


CJOUT C. 

? .Cvu 

C.O 

LB /SEC 


** ******* **************** 

********** ***** 


DJCT 


- 

' ‘ * ■ ' — " ■ 


1 NO VEw T TOT 

P TOT 

P STAT AREA 

GAN 

.200 434 . 2100. 

16151. 

1573 A 

. .1990 

1.3130 

JTIPKAX STRESS 

DEN 

W /AREA 

TR H/T 


1205.2 1 i 4 2 5 . 5 

• 2 56 

.547 

5.000 

748 






> 


TURBINE 

4 MECHANICAL 

DESIGN 




H/T 

N STG 

LOADING 

AREA 




• 74b 

3 . C f 0 

• 4 2 0 

•T9 9 

- 



UT 

KTIP 

RHUB 

del h 

RPM 

HAXRPM 

TCRQ 

1 2 3 5.2 

4.6 

3.4 

207.3 

30352.6 

30352.6 

7429 


Figure 44. (continued) 
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Figure 44. ( continued, ) 
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1AX CONDITIONS OLCUR hT 


H \ 

r. • a j 


ALT 

PTCT 

ttot c. 

********* 

NOZZLE 5 

*^1GHT= 19.73 LENGTH 


*> • U L. J 

******** 


14.9ics TR WT = 


0 • C u 


A ************ 
* * 

* SHAF 6 * 

* * 
************2 


iax torjui condition 

*********************** 
TORQ JE 

• 9 

*********************** 

SHAFT f» 

03 31 LENlG 

1.&5 l«4-.i 13.32 


ON 

1.27 


UT 

2 .37 


T3TAL INERTIA 0 r THIS SPOOL IS 13A8. 


************* 
* * 

* ACCS WT * 




ACCS T = *» . ," 3 


Figure 44 . (continued) 
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Tat** 9 Error 


1. 

" ** f uv 

Compressor, I, pressure ratio is too high" - more than 20 

compressor stages calculated. First stage maximum pressure 
ratio too small. 

2. 

Compressor, I. stage and blade parameters, meaningless" - 
stage inlet Mach number less than or equal to zero, or hub 
radius of compressor equals zero. 

3. 

"Duct is not converging - error only called for rotating 
splitter fan component. Inlet or exit Mach numbers of fan 
may be input incorrectly. 

4. 

"Frror in shaft" - iteration for shaft diameter not converging. 
Check shaft inputs. 

5. 

"Turbine, I, work or radius too high. RC = , X.XX" - more 

than 9 turbine stages calculated - turbine loading parameter 
too small or control radius improperly input. 

6. 

Turbine, I, stage and blade parameters meaningless" - Mach 
number or hub radius less than or equal to zero 
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?.6 Program Structure 


I 


The execution flow was designed to minimize the interaction of the basic 
NNEP program and the weight estimation routines. The only data flowing 
between the two is via the common blocks SNGL and DBL and the variables 
IWT and IPLT. The subroutine THERM is used to obtain thermodynamic 

properties of the fluid. An assumption is made that the thermodynamic 
properties are established at each station prior to calling WTEST 
subroutine, which is flow charted in Figure 45. 

WTEST calls the component routines. These component routines are 
independent of each other, and some use the same lower level routines as 
some other components. After all weicihts and dimensions have been 

estimated ENGPLT is called to make the printer plot. For a description 
of subroutine connectivity, see Figure 46. 

The code is in FORTRAN IV and has been checked out on IBM 370/168. The 
code is single precision except for the values in the Navy-NASA Engine 
Program (NNEP). The code was designed to minimize conversion 

requirements to other machines. No subroutines are required beyond those 

in the IBM FORTRAN IV manual, and there is no character manipulation, 

only full word tests are used when testing BCD input. 

The NNEP/WATE-1 code requires 75B49^ core (482120 10 bytes) to run 
without buffers. The execution of a design point followed by a weight 

estimation and printer plot is 4 seconds CPU time, using an existinq load 
module (all data reference 370/168). 

The fo! low inn variables in NNEP common blocks may be referenced by a 
component, weiqht estimating routine depending on the component type: 
DATOUT, WTF , TOPRES, TOTEMP, FAR, CORFLO, JCONF , JTYPE, NCOMP, NOSTAT, 
NF INIS. In no case is any value changed by the weiqht estimation code. 
Each call to WTEST routine will cause a N/\MFLIST read of "W" data. This 
is the one and onlv read in weiqht estimation code. Based on the 

information in NCOMP and JTYPE the proper component routine is called 
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Figure 45. Functional Flow Chart of WTEST 
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with the competent number (I) as an argument. Each component is expected 
to fill WATE (n, ALF.NG fl), TLENG (I), RO H,!), RO F2,Il, RI (1,11, RI 
|?,n. Rota*:inq components also fill RPMT (II. The shaft component fills 
OSHAFT fNl where N is the shaft count from the inside out. 




The array CONVER in common CONVER are conversion factors to convert English units to SI 


units: 


ARRAY# 

VALUE 

UNITS 

ENGLISH 

to 

UNITS 

SI 

1 

2.54 

inch 


cin 


.3048 

feet 


meter 

3 

453b 



Kg 

4 

.0020 

ft- 


meter 3 

5 

.02768 

ib m /iu 3 


K g /cm 3 

b 

.680475 

Ibj/in- 


Newton/cm 3 

1 

4.882 



Kg/m- 

8 

.555 

°R 


°K 

9 

1.05435 

BTU/sec 


K WATTS 

10 

.07457 

HP 


K WATTS 

1 1 

47.88 

Ibf/ft- 


Newton/m ~ 



— 0 CONCLUSIONS AND RECOMMENDATIONS 


The WATE-2 method can provide a reasonable level of accuracy in 
predicting weight and dimensions of large and small gas turbine engines. 
Flexibility of modeling the engine cycle and gas-path allows virtually 
any conventional or non-conventional gas turbine engine to be analyzed. 
Normally, all components can be sized from the stage by stage calculation 
of state conditions that are accomplished internally. 

A number of new capabilities have been added to the original program that 
will also make it more useful: small gas turbines, rotating interia, 

airflow scaling, center of gravity, engine design diagnostic aids, and 
flight envelope maximization of component size. Accuracy of the program 
has also been improved with the addition of a new preliminary design 
procedure for discs. The disc weight correlations that were previously 
developed did not include the effects of blade aspect ratio or disc 
stress; it was onlv sensitive to blade pull stress and disc radius. Both 
methods are available as options, however it is recommended that the new 

disc design procedure be used for all types of engines rather than the 

disc weight correlation procedure. 

Further improvements to the program are also recommended. An engine 
desion procedure similar to that described in Reference 8 could be built 
into the program to eliminate external iterations of the engine design. 
Exit Mach number in turbines should he used to establish blade sizes 
rather than the entry Mach number as now used. Since there is usually 

less variation in the exit angle and exit Mach number than there is at 

the entry of turbines, this would provide a better blade sizing 
criteria. The program could also be changed to a constant-mean radius 
design procedure rather than the constant-area area definition that is 
currently used. The former method is more cormonlv used in the industry. 
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